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Abstract 

These lectures provide an overview of Quantum Chromodynamics 
(QCD), the SU(3)c gauge theory of the strong interactions. Af- 
ter briefly reviewing the empirical considerations which lead to the 
introduction of colour, the QCD Lagrangian is discussed. The run- 
ning of the strong coupling and the associated property of Asymp- 
totic Freedom are analyzed. Some selected experimental tests and 
the present knowledge of a s are summarized. A short description 
of the QCD flavour symmetries and the dynamical breaking of 
chiral symmetry is also given. A more detailed discussion can be 
found in standard textbooks and recent reviews 



1. QUARKS AND COLOUR 

A fast look into the Particle Data Tables [[HJ reveals the richness and variety of 
the hadronic spectrum. The large number of known mesonic and baryonic states clearly 
signals the existence of a deeper level of elementary constituents of matter: quarks 
In fact, the messy hadronic world can be easily understood in terms of a few constituent 
spin- 1 quark flavours: 
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Assuming that mesons are M = qq states, while baryons have three quark constituents, 
B = qqq, one can nicely classify the entire hadronic spectrum: 
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*) Lectures given at the 1994 European School of High-Energy Physics (Sorrento, Italy, 29 
August - 11 September, 1994) 
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There is a one-to-one correspondence between the observed hadrons and the states 
predicted by this simple classification; thus, the Quark Model appears to be a very useful 
Periodic Table of Hadrons. However, the quark picture faces a problem concerning the 
Fermi-Dirac statistics of the constituents. Since the fundamental state of a composite 
system is expected to have L = 0, the A ++ baryon (J = |) corresponds to u^u^v) , with 
the three quark-spins aligned into the same direction (53 = +|) and all relative angular 
momenta equal to zero. The wave function is symmetric and, therefore, the A ++ state 
obeys the wrong statistics. 



The problem can be solved assuming |12| the existence of a new quantum number, 
colour, such that each species of quark may have No = 3 different colours: q a , a = 1, 2, 3 
(red, yellow, violet). Then one can reinterpret the A ++ as the antisymmetric state 
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(notice that at least 3 colours are needed for making an antisymmetric state). In this 
picture, baryons and mesons are described by the colour-singlet combinations 
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In order to avoid the existence of non-observed extra states with non-zero colour, one needs 
to further postulate that all asymptotic states are colourless, i.e. singlets under rotations in 
colour space. This assumption is known as the confinement hypothesis, because it implies 
the non-observability of free quarks: since quarks carry colour they are confined within 
colour-singlet bound states. 

The quark picture is not only a nice mathematical scheme to classify the hadronic 
world. We have strong experimental evidence of the existence of quarks. Fig. [I] shows a 
typical Z — > hadrons event, obtained at LEP. Although there are many hadrons in the 
final state, they appear to be collimated in 2 jets of particles, as expected from a two-body 
decay Z — > qq, where the qq pair has later hadronized. 




Figure 1: Two-jet event from the hadronic decay of a Z boson (DELPHI). 
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Figure 2: Feynman diagram for e + e — > hadrons. 



1.1 Evidence of colour 

A direct test of the colour quantum number can be obtained from the ratio 

a(e + e~ — > hadrons) 

R e +e- = / 7 _ 7 _n • (1-3) 

cr(e+e — > n + n ) 

The hadronic production occurs through e + e~ — > 7*, Z* — >• qq — > hadrons. Since quarks 
are assumed to be confined, the probability to hadronize is just one; therefore, the sum 
over all possible quarks in the final state will give the total inclusive cross-section into 
hadrons. At energies well below the Z peak, the cross-section is dominated by the 7- 
exchange amplitude; the ratio R e + e - is then given by the sum of the quark electric charges 
squared: 

(Nf = 3 : u,d,s) 

Re + e-^N c J2Q'j = {fN c = f, (N f = A : u,d,s,c) . (1.4) 

(Nf = 5 : u,d,s, c, b) 
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Figure 3: Measurements of R e + e - 
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The two continuous curves are QCD fits. 



Figure 4: r-decay diagram. 



The measured ratio is shown in Fig. [| Although the simple formula (|1.4j ) cannot 



explain the complicated structure around the different quark thresholds, it gives the right 
average value of the cross-section (away from the thresholds), provided that Nc is taken 
to be three. The agreement is better at larger energies. Notice that strong interactions 
have not been taken into account; only the confinement hypothesis has been used. 

The hadronic decay of the r lepton provides additional evidence for Nc = 3. The 
decay proceeds through the ly-emission diagram shown in Fig. f|. Since the W coupling 
to the charged current is of universal strength, there are (2 + Nc) equal contributions (if 
final masses and strong interactions are neglected) to the r-decay width. Two of them 
correspond to the leptonic decay modes r~ — > z/ T e~z/ e and t~ — ► Urfi^u^, while the other 
Nc are associated with the possible colours of the quark-ant iquark pair in the t~ — > u T dgu 
decay mode (de = cos 9cd+ sin 9c s). Hence, the branching ratios for the different channels 
are expected to be approximately: 

B„, = Bt(t" - v T n>,) « — L- = i = 20% , (1.5) 
T(t -> u T e~u e ) 



which should be compared with the experimental averages [10]: 



B T ^ e = (18.01 ±0.18)%, B T ^ = (17.65 ±0.24)%, (1.7) 

R T = (1 - B T _+ e - B T ^)/B T ^ e = 3.56 ± 0.04 . (1.8) 

The agreement is fairly good. Taking Nc = 3, the naive predictions only deviate from the 
measured values by about 20%. Many other observables, such as the partial widths of the 
Z and bosons, can be analyzed in a similar way to conclude that Nc = 3. 

A particularly strong test is obtained from the tt° — > 77 decay, which occurs through 
the triangular quark loops in Fig. |5|. The crossed vertex denotes the axial current A 3 ^ = 
(uj^u - dj^d). One gets: 

(39 , sS-7.7».v. M 

where the 7r° coupling to A 3 , f n = 92.4 MeV, is known from the 7r~ — > /i _ z/ M decay rate 
(assuming isospin symmetry). The agreement with the measured value, T = 7.7 ± 0.6 
eV [0, is remarkable. With Nc = 1, the prediction would have failed by a factor of 9. 



The nice thing about this decay is that it is associated with an anomaly: a global flavour 
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Figure 5: Triangular quark loops generating the decay tt° — > 77. 

symmetry which is broken by quantum effects (the triangular loops). One can then proof 
that the decay amplitude ( |L9p does not get corrected by strong interactions [T3| . 

Anomalies provide another compelling theoretical reason to adopt Nc = 3. The 
gauge symmetries of the Standard Model of electroweak interactions have also anomalies 
associated with triangular fermion loops (diagrams of the type shown in Fig. |5|, but with 
arbitrary gauge bosons -W ± , Z, 7- in the external legs and Standard Model fermions in 
the internal lines). These gauge anomalies are deathly because they destroy the renor- 
malizability of the theory. Fortunately, the sum of all possible triangular loops cancels 
if Nq = 3. Thus, with three colours, anomalies are absent and the Standard Model is 
well-defined. 

1.2 Asymptotic Freedom 

The structure of the proton can be probed through the scattering e~p — > e~p. The 
cross-section is given by 



da TroWf jl^(g 2 )| 2 + ^|G M (Q 2 )| 2 Q 2 Q 
' ) ' +77T7^M(Q 2 )| 2 tan 2 - 



dQ 2 AE 2 sin 4 ~EE' ] 1 + ^ 2M p 



;i.io) 



where E and E' are the energies of the incident and scattered electrons, respectively, in 
the proton rest-frame, 9 the scattering angle, M p the proton mass and 

Q 2 = -q 2 = AEE'sm 2 ^, (1.11) 

with g M = (k e — k' e ) )l the momentum transfer through the intermediate photon propagator. 
Gg and Gm 

are the electric and magnetic form factors, respectively, describing the 
proton electromagnetic structure; they would be equal to one for a pointlike spin-^ target. 
Experimentally they are known to be very well approximated by the dipole form 

G M (Q 2 )/^ P « G E m + ^-^) , (1-12) 

where \x v = 2.79 is the proton magnetic moment (in proton Bohr magneton units). Thus, 
the proton is actually an extended object with a size of the order of 1 fm. At very low 
energies (Q 2 « 1 GeV 2 ), the photon probe is unable to get information on the proton 
structure, Gm,e(Q 2 ) ~ Gm,e(0) = 1, and the proton behaves as a pointlike particle. At 
higher energies, the photon is sensitive to shorter distances; the proton finite size gives 
then rise to form factors, which suppress the elastic cross-section at large Q 2 , i.e. at large 
angles. 
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Figure 6: Inelastic e p — > e X scattering. 



One can try to further resolve the proton structure, by increasing the incident 
energy. The inelastic scattering e~p — ► e~X becomes then the dominant process. Making 
an inclusive sum over all hadrons produced, one has an additional kinematical variable 
corresponding to the final hadronic mass, W 2 = P\. The scattering is usually described 
in terms of Q 2 and 

M p 2M p ' v ; 

where P M is the proton cuadrimomentum; v is the energy transfer in the proton rest-frame. 
In the one-photon approximation, the unpolarized differential cross-section is given by 

da 7ra; 2 cos 2 ^ f . . ,9. <>0l 

V 2 (Q 2 ,z/) + 2iy 1 (Q 2 ,z/) tan 2 - . (1.14) 
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dQ 2 dv 4£ 2 sin 4 §££ 
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The proton structure is then characterized by two measurable structure functions. For a 
pointlike proton, the elastic scattering ( |1.10| ) corresponds to 

w ^ = •*<«"■">='("-§;)■ (hl5) 

At low Q 2 , the experimental data reveals prominent resonances; but this resonance 
structure quickly dies out as Q 2 increases. A much softer but sizeable continuum contri- 
bution persists at large Q 2 , suggesting the existence of pointlike objects inside the proton. 

To get an idea of the possible behaviour of the structure functions, one can make 
a very rough model of the proton, assuming that it consist of some number of pointlike 
spin- 1 constituents (the so-called partons), each one carrying a given fraction ^ of the 
proton momenta, i.e. p^ = ^iP^. That means that we are neglecting^ the transverse 
parton momenta, and rrii = £M P . The interaction of the photon-probe with the parton i 
generates a contribution to the structure functions given by: 

W?{Q\v) = ^s(u-gr)=^5&-x), (1.16) 




Wf(QV) = e 2 s[u-^-)=e 2 ^-6U t -x), (1.17) 



These approximations can be made more precise going to the infinite momentum frame of 
the proton, where the transverse motion is negligible compared with the large longitudinal 
boost of the partons. 
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where e; is the parton electric charge and 

Q 2 



x = 



Q 1 



2M p u Q 2 + W 2 - Ml 



:i-18) 



Thus, the parton structure functions only depend on the ratio x, which, moreover, fixes 
the momentum fractions We can go further, and assume that in the limit Q 2 — > oo, 
v — > oo, but keeping x fixed, the proton structure functions can be estimated from an 
incoherent sum of the parton ones (neglecting any strong interactions among the partons). 
Denoting /;(£«) the probability that the parton % has momentum fraction £j, one then has: 



X 
V 



(1.19) 
(1.20) 



This simple parton description implies then the so-called Bjorken scaling [f[4j| : the proton 
structure functions only depend on the kinematical variable x. Moreover, one gets the 
Callan-Gross relation |15| 

F 2 {x) = 2xF 1 {x) , (1.21) 

which is a consequence of our assumption of spin-| partons. It is easy to check that spin-0 
partons would have lead to Fi(x) = 0. 




0.2 0.4 0.6 0.8 




Figure 8: The ratio 2xFi/F 2 versus x, for 
different Q 2 values (1.5 GeV 2 < Q 2 < 16 
GeV 2 ) H (taken from Ref. 0). 



Figure 7: Experimental data on vW 2 as 
function of x, for different values of Q 2 [fL6| 
(taken from Ref. (l[|). 

The measured values of vW 2 (Q 2 , v) are shown in Fig. [7| as function of x, for many 
different values of Q 2 between 2 and 18 GeV 2 ; the concentration of data points along a 
curve indicates that Bjorken scaling is correct, to a quite good approximation. Fig. |8] shows 
that the Callan-Gross relation is also reasonably well satisfied by the data, supporting 
the spin- 1 assignment for the partons. 
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The surprising thing of this successful predictions is that we have assumed the 
existence of free independent pointlike partons inside the proton, in spite of the fact that 
quarks are supossed to be confined by very strong colour forces. Bjorken scaling suggests 
that the strong interactions must have the property of asymptotic freedom: they should 
become weaker at short distances, so that quarks behave as free particles for Q 2 — > oo. 
This also agrees with the empirical observation in Fig. [3|, that the free-quark description 
of the ratio R e + e - works better at higher energies. 

Thus, the interaction between a qq pair looks like some kind of rubber band. If we 
try to separate the quark form the antiquark the force joining them increases. At some 
point, the energy on the elastic band is bigger than 2m q >, so that it becomes energetically 
favourable to create an additional q'q' pair; then the band breaks down into two mesonic 
systems, qq' and q'q, each one with its corresponding half-band joining the quark pair. 
Increasing more and more the energy, we can only produce more and more mesons, but 
quarks remain always confined within colour-singlet bound states. Conversely, if one tries 
to approximate two quark constituents into a very short-distance region, the elastic band 
loses the energy and becomes very soft; quarks behave then as free particles. 



1.3 Why SU{3)7 

Flavour-changing transitions have a much weaker strength than processes mediated 
by the strong force. The quark-flavour quantum number is associated with the electroweak 
interactions, while strong forces appear to be flavour-conserving and flavour- independent. 
On the other side, the carriers of the electroweak interaction (7, Z, W^) do not couple 
to the quark colour. Thus, it seems natural to take colour as the charge associated with 



the strong forces and try to build a quantum field theory based on it [18|]. The empirical 
evidence described so far puts a series of requirements that the fundamental theory of 
colour interactions should satisfy: 

1. Colour is an exact symmetry Gq (hadrons do not show colour multiplicity). 

2. N c = 3. Thus, quarks belong to the triplet representation 3 of Gc- 

3. Quarks and antiquarks are different states. Therefore, 3* 7^ 3, i.e. the triplet repre- 
sentation has to be complex. 

4. Confinement hypothesis: hadronic states are colour singlets. 

5. Asymptotic freedom. 

Among all compact simple Lie groups there are only four having 3-dimensional 
irreducible representations; moreover, three of them are isomorphic to each other. Thus, we 
have only two choices: SU(3) or 50(3) ~ SU(2) ~ Sp(l). Since the triplet representation 
of S0(3) is real, only the symmetry group SU(3) survives the conditions 1, 2 and 3. The 
well-known SU(3) decomposition of the products of 3 and 3* representations, 



qq 

qqq 

qq 

qqqq 



3® 3* = 1©8, 

3©3©3 = 1 ©8©8© 10, 

3©3 = 3*©6, 

3©3©3©3 = 3©3©3©6*©15©15©15©15:, (1.22) 



guarantees that there are colour-singlet configurations corresponding to meson (qq) and 
baryon (qqq) states, as required by the confinement hypothesis. Other exotic combinations 
such as diquarks (qq, qq) or four-quark states (qqqq, qqqq) do not satisfy this requirement. 
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Clearly, the theory of colour interactions should be based on the SU(3)c group. It 
remains to be seen whether such a theory is able to explain confinement and asymptotic 
freedom as natural dynamical consequences of the colour forces. 

2. GAUGE SYMMETRY: QED 

Let us consider the Lagrangian describing a free Dirac fermion: 

C = i^{x)Ydp^{x) - m^{x)^{x). (2.1) 
Cq is invariant under global U(l) transformations 

V'(x) = exp{iQ9}$(x), (2.2) 

where Q9 is an arbitrary real constant. The phase of \fr(x) is then a pure convention- 
dependent quantity without physical meaning. However, the free Lagrangian is no-longer 
invariant if one allows the phase transformation to depend on the space-time coordinate, 
i.e. under local phase redefinitions 9 = 9(x), because 

d^(x) ^ exp {iQ6} (dp + iQdp9) *(x) . (2.3) 

Thus, once an observer situated at the point Xq has adopted a given phase-convention, 
the same convention must be taken at all space-time points. This looks very unnatural. 

The "Gauge Principle" is the requirement that the U(l) phase invariance should 
hold locally. This is only possible if one adds some additional piece to the Lagrangian, 
transforming in such a way as to cancel the dp9 term in Eq. ( |2.3| ). The needed modification 
is completely fixed by the transformation ( |2.3| ): one introduces a new spin-1 (since dp9 
has a Lorentz index) field Ap(x), transforming as 

Ap(x) ^ A'^x) ee Ap(x) + - e 8,9, (2.4) 

and defines the covariant derivative 

D^(x) = [8^ - ieQA^(x)] , (2.5) 
which has the required property of transforming like the field itself: 

D^(x) ^ (DpV)'(x) = exp{iQ9}D^(x). (2.6) 

The Lagrangian 

C = i^(x)YD^{x) - m*(z)*(x) = C + eQA lt (x)V(x)'f*(x) (2.7) 

is then invariant under local U(l) transformations. 

The gauge principle has generated an interaction between the Dirac spinor and 
the gauge field Ap, which is nothing else than the familiar QED vertex. Note that the 
corresponding electromagnetic charge eQ is completely arbitrary. If one wants A^ to be 
a true propagating field, one needs to add a gauge-invariant kinetic term 

^Kin = —-rFpfF^, (2.8) 
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where F^ u = d^A u — d u A^ is the usual electromagnetic field strength. A possible mass 
term for the gauge field, ^m 2 A^A )1 , is forbidden because it would violate gauge invariance; 
therefore, the photon field is predicted to be massless. The total Lagrangian in ( |2.7| ) and 
(|2.8| ) gives rise to the well-known Maxwell equations. 

From a simple gauge- symmetry requirement, we have deduced the right QED La- 
grangian, which leads to a very successful quantum field theory. Remember that QED 
predictions have been tested to a very high accuracy, as exemplified by the electron and 



muon anomalous magnetic moments [a/ = (gi — 2)/2, where /i; = g\ (e%/2mi)\ [19]: 

(115 965 214.0 ± 2.8) x KT 11 (Theory) 
(115 965 219.3 ± 1.0) x 10~ n (Experiment)' 



, (1 165 919.2 ± 1.9) x 10~ 9 (Theory) 
a„ = < „ . 2.10) 

1 (1 165 923.0 ± 8.4) x 10~ 9 (Experiment) 



3. THE QCD LAGRANGIAN 

Let us denote qj a quark field of colour a and flavour /. To simplify the equa- 
tions, let us adopt a vector notation in colour space: qj = column(gj, q"j, q'j) . The free 
Lagrangian 

£ o = Vf (^ d v - m f)Qf (3- 1 ) 
/ 

is invariant under arbitrary global SU(3)c transformations in colour space, 

Qf — > (if)' = U a p q P } , UU ] = U ] U = 1 , det U = 1 . (3.2) 

The SU(3)c matrices can be written in the form 

U = exp |-Z3> s y0 a j , (3.3) 

where A a (a = 1, 2, . . . , 8) denote the generators of the fundamental representation of 
the SU(3)c algebra, and 6 a are arbitrary parameters. The matrices A a are traceless and 
satisfy the commutation relations 



X\\ b 



2if abc X c , (3.4) 



with f abc the SU(3)c structure constants, which are real and totally antisymmetric. Some 
useful properties of SU(3) matrices are collected in Appendix A. 

As in the QED case, we can now require the Lagrangian to be also invariant under 
local SU{3)c transformations, 6 a = 9 a (x). To satisfy this requirement, we need to change 
the quark derivatives by covariant objects. Since we have now 8 independent gauge pa- 
rameters, 8 different gauge bosons G^(x), the so-called gluons, are needed: 



D»q f = 



q } = [df-igsG^x)] q f . (3.5) 



Notice that we have introduced the compact matrix notation 



[G»{x)U = ( y J G»(x) . (3.6) 
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We want D^qj to transform in exactly the same way as the colour- vector qj; this fixes 
the transformation properties of the gauge fields: 

D n — , ( D ny = UD »u\. — > f G ny = jj jj\ — 1 _ (^m _ ( 3 _7) 

9s 

Under an infinitesimal SU(3)c transformation, 

qj — W = <ff-ig.(jJ a 89 aq ;, 

G£ — (GZy = GZ-d»(56 a )+g s f abc 56 b GZ. (3.8) 

The gauge transformation of the gluon fields is more complicated that the one obtained 
in QED for the photon. The non-commutativity of the SU(3)c matrices gives rise to an 
additional term involving the gluon fields themselves. For constant 58 a , the transforma- 
tion rule for the gauge fields is expressed in terms of the structure constants f abc only; 
thus, the gluon fields belong to the adjoint representation of the colour group (see Ap- 
pendix A). Note also that there is a unique SU (3)c coupling g s . In QED it was possible to 
assign arbitrary electromagnetic charges to the different fermions. Since the commutation 
relation ( |3.4| ) is non-linear, this freedom does not exist for SU{3)c- 

To build a gauge-invariant kinetic term for the gluon fields, we introduce the cor- 
responding field strengths: 

G^ix) = — [D", D v ] = 3 fl G u - d u G» - ig s [G", G v \ = — G% v (x) , 
9s 2 

G?(x) = d'Gl-d-GZ + gJ^GZG". (3.9) 
Under a gauge transformation, 

— >{G^)' = UG^U\ (3.10) 

and the colour trace r Yx(G y ' v G ^v) = \G^ u G a ^ u remains invariant. 

Taking the proper normalization for the gluon kinetic term, we finally have the 
SU(3)c invariant QCD Lagrangian: 

£qcd = -\ G^G% + J2lf (rfDn - m f ) q f . (3.11) 

It is worth while to decompose the Lagrangian into its different pieces: 

£qcd = -\{d»Gl-d»G»){d»Gl-d v Gl)+Y,$ W d »~™f) 

f 

+ 9.<%E%k(t) q f (3 - 12) 

_ | fabc {d , G u _ d u G ,j _ 9j fabcfade G £ G:G d G e 

The first line contains the correct kinetic terms for the different fields, which give rise to 
the corresponding propagators. The colour interaction between quarks and gluons is given 
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by the second line; it involves the SU(3)c matrices A a . Finally, owing to the non-abelian 
character of the colour group, the G^G a term generates the cubic and quartic gluon 
self-interactions shown in the last line; the strength of these interactions is given by the 
same coupling g s which appears in the fermionic piece of the Lagrangian. 

In spite of the rich physics contained in it, the Lagrangian (ETO looks very simple, 
because of its colour-symmetry properties. All interactions are given in terms of a single 
universal coupling g s , which is called the strong coupling constant. The existence of self- 
interactions among the gauge fields is a new feature that was not present in the QED 
case; it seems then reasonable to expect that these gauge self-interactions could explain 
properties like asymptotic freedom and confinement, which do not appear in QED. 




Figure 9: Three-jet event from the hadronic decay of a Z boson (DELPHI). 

Without any detailed calculation, one can already extract qualitative physical con- 
sequences from £qcd- Quarks can emit gluons. At lowest-order in g s , the dominant process 
will be the emission of a single gauge boson. Thus, the hadronic decay of the Z should 
result in some Z — > qqG events, in addition to the dominant Z — > qq decays discussed 
in Section |l] Fig. |^ clearly shows that 3-jet events, with the required kinematics, indeed 
appear in the LEP data. Similar events show up in e + e~ annihilation into hadrons, away 
from the Z peak. 

In order to properly quantize the QCD Lagrangian, one needs to add to £qcd the 
so-called Gauge-fixing and Faddeev-Popov terms. Since this is a rather technical issue, its 
discussion is relegated to Appendix B. 



4. QUANTUM LOOPS 

The QCD Lagrangian is rather economic in the sense that it involves a single 
coupling g s . Thus, all strong-interacting phenomena should be described in terms of just 
one parameter. At lowest order in g s (tree-level), it is straightforward to compute all 
kind of scattering amplitudes involving quarks and gluons: qq — ► GG, qq — > qq, Gq — ► 
Gq, . . . Unfortunately, this exercise by itself does not help very much to understand the 
physical hadronic world. First, we see hadrons instead of quarks and gluons. Second, 
we have learnt from experiment that the strength of the strong forces changes with the 
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energy scale: the interaction is very strong (confining) at low energies, but quarks behave 
as nearly free particles at high energies. Obviously, we cannot understand both energy 
regimes with a single constant g s , which is the same everywhere. Moreover, if we neglect 
the quark masses, the QCD Lagrangian does not contain any energy scale; thus, there is 
no way to decide when the energy of a given process is large or small, because we do not 
have any reference scale to compare with. 

If QCD is the right theory of the strong interactions, it should provide some dy- 
namical scale through quantum effects. 



4.1 Regularization of loop integrals 

The computation of perturbative corrections to the tree-level results involves diver- 
gent loop integrals. It is then necessary to find a way of getting finite results with physical 
meaning from a priori meaningless divergent quantities. 




a n x ^__^ / v b 
k-q 

Figure 10: Gluon self-energy diagram. 



Let us consider the self-energy gluon loop in Fig. 10. The corresponding contribu- 



tion in momentum space can be easily obtained, using standard Feynman rules techniques: 

Tj^f \ 2cm/ d*k Tr[7^rg ~ 4)} (AU 

The result is proportional to g 2 , because there are two qqG vertices, and there is a trivial 
SU(3)c factor, Tp = ~, coming from the colour trace |Tr(A a A b ) = 5 ab T F . 

The problem appears in the momentum integration, which is clearly divergent [~ 
J d 4 k(l/k 2 ) = oo]. We can define the momentum integral in many different (and arbitrary) 
ways. For instance, we could introduce a cut-off M, such that only momentum scales 
smaller than M are integrated; obviously, the resulting integral would be an increasing 
function of M. Instead, it has become conventional to define the loop integrals through 
dimensional regularization: the calculation is performed in D = 4 + 2e dimensions. For 
e 7^ the resulting integral is well-defined: 

d D k k^k-qf _ -i f-q'Y^ ( x 5 \ f gV* | A ^ (42) 



(2tt) d k 2 (k-q) 2 6(4vr) 2 \ 4tt J v ; V 3/[2(l + e 



The ultraviolet divergence of the loop appears at e = 0, through the pole of the Gamma 
function, 

r(-e) = -I- 7E + 0( e 2 ), (4.3) 
where jp = 0.577215 ... is the Euler constant. 
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Although the integral ( |4.2| ) looks somewhat funny, dimensional regularization has 
many advantages because does not spoil the gauge symmetry of QCD and, therefore, 
simplifies a lot the calculations. One could argue that a cut-off procedure would be more 
physical, since the parameter M could be related to some unknown additional physics at 
very short distances. However, within the QCD framework, both prescriptions are equally 
meaningless. One just introduces a regularizing parameter, such that the integral is well- 
defined and the divergence is recovered in some limit (M — > oo or e — > 0). 

Since the momentum-transfer q 2 has dimensions, it turns out to be convenient to 
introduce and arbitrary energy scale /i and write 



-T 



:2c 



-r 



:2c 



Ann 2 



-e) = -/x 



+ lE - ln47r + In (-g 2 //x 2 ) + 0(e) 



(4.4) 

Obviously, this expression does not depend on /x; but written in this form one has a 
dimensionless quantity {—q 2 /fJ> 2 ) inside the logarithm. 

The contribution of the loop diagram in Fig. 1C can finally be written as 



11 ab 



U(q 2 



Sab 
4 



-g 2 <T 



9sV 



q»q u 



i% 2 



47T 



+ lE - ln4vr + In (-g 2 //x 2 ) - \ + 0(e)} . (4.5) 



Owing to the ultraviolet divergence, Eq. ( |4.5|) does not determine the wanted self- 
energy contribution. Nevertheless, it does show how this effect changes with the energy 
scale. If one could fix the value of IT(g 2 ) at some reference momentum transfer q$, the 
result would be known at any other scale: 

2 



n(^) = n(^) - -t f 



(4.6) 



We can split the self-energy contribution into a meaningless divergent piece and a 
finite term, which includes the q 2 dependence, 



U(q 2 ) = An e (/i 2 )+n fl (g 2 //i 2 



(4.7) 



This separation is of course ambiguous, because the finite g 2 -independent contributions 
can be splitted in many different ways. A given choice defines a scheme: 



ail(/x 2 ) 



_2> 9i..2e 
3tt 47r^ 

3?r 47T ' e 

3tt 4tt^ 

07T 4"7T 



7s - ln(47r) 



+ 7s - ln(4vr) 



.Ik Si 

3n 4"7T 

.If si 

37Y 47T 



-<z 2 /V) 



+ - ln(47r) 



(/x-scheme), 

(MS-scheme), (4.8) 
(MS-scheme), 

(/i-scheme), 

(MS-scheme), (4.9) 
(MS-scheme) . 



In the /i-scheme, one uses the value of II(— /x 2 ) to define the divergent part. MS and 
MS stand for minimal subtraction [[20] and modified minimal subtraction schemes [21 



in the MS case, one subtracts only the divergent 1/e term, while the MS scheme puts 
also the annoying 7^ — ln(47r) factor into the divergent part. Notice that the logarithmic 
^-dependence is always the same. 
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Figure 11: Photon self-energy contribution to e e scattering. 

4.2 Renormalization: QED 

A Quantum Field Theory is called renormalizable if all ultraviolet divergences can 
be reabsorbed through a redefinition of the original fields and couplings. 

Let us consider the electromagnetic interaction between two electrons. At one loop, 
the QED photon self-energy contribution is just given by Eq. [4.5| , with the changes Tp — > 1 
and g s — > e. The corresponding scattering amplitude takes the form 

Ttf) „ -J^t {l_n(g 2 ) + ...} , (4.10) 

where J M denotes the electromagnetic fermion current. 

At lowest order, T(q 2 ) ~ ctjq 2 with a = e 2 /(4n). The divergent correction gener- 
ated by quantum loops can be reabsorbed into a redefinition of the coupling: 

^ {1 - AITV) - lW//i 2 )} - ^ {1 - IW//^ 2 )} , (4.11) 

a R (fi 2 ) = a (l + ^/i 2e [- + C schem el +...) , a = ^ , (4.12) 

I 3n Le J J 47r 

where eo denotes the bare coupling appearing in the QED Lagrangian; this bare quantity is, 

however, not directly observable. Making the redefinition ( |4.11[ ), the scattering amplitude 



is finite and gives rise to a definite prediction for the cross-section, which can be compared 
with experiment; thus, one actually measures the renormalized coupling a«. 

The redefinition ( [4.1 1|) is meaningful, provided that it can be done in a self- 
consistent way: all ultraviolet divergent contributions to all possible scattering processes 
should be eliminated through the same redefinition of the coupling (and the fields). The 
nice thing of gauge theories, such as QED or QCD, is that the underlying gauge symmetry 
guarantees the renormalizability of the quantum field theory. 

The renormalized coupling a R (n 2 ) depends on the arbitrary scale fi and on the 
chosen renormalization scheme [the constant C sc heme denotes the different finite terms 
in Eq. (|4.8|) ]. Quantum loops have introduced a scale dependence in a quite subtle way. 
Both a R ([i 2 ) and the renormalized self-energy correction U R (q 2 / fi 2 ) depend on fi, but the 
physical scattering amplitude T(q 2 ) is of course //-independent: (Q 2 = —q 2 ) 



T(q 2 ) 4ttJ%^P 1 + 



q 2 I 3tt 



,2 s 



^ I o J ^scheme 



+ 



47T J% —gT~ I 1 + 3n C scheme + •"> ■ (4-13) 
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The quantity ot(Q 2 ) = an(Q 2 ) is called the QED running coupling. The ordinary 
fine structure constant a = 1/137 is defined through the classical Thomson formula; 
therefore, it corresponds to a very low scale Q 2 = —m 2 . Clearly, the value of a relevant 
for LEP experiments is not the same [o;(M|)j^g = 1/129]. The scale dependence of a{Q 2 ) 
is regulated by the so-called /3-function: 



H— = a /3(a); /9(a) = fa- + (3 2 - 

Ct/X 71 V 7T / 



+ ■ 



(4.14) 



At the one- loop level, the (3-i unction reduces to the first coefficient, which is fixed by 
Eq. (|4121) : 

; QED _ 2 

3 ' 



ft 



(4.15) 



The first-order differential equation ( |4.14|) can then be easily solved, with the result: 



a(Q 2 



a(«g) 



2tt 



In (Q 2 /Qo) 



(4.16) 



Since /3i > 0, the QED running coupling increases with the energy scale: a(Q 2 ) > q(Qq) 
if Q 2 > Ql; i.e. the electromagnetic charge decreases at large distances. This can be 
intuitively understood as the screening effect due to the virtual e + e~ pairs generated, 
through quantum effects, around the electron charge. The physical QED vacuum behaves 
as a polarized dielectric medium. 




Figure 12: Electromagnetic charge screening in a dipolar medium. 



Notice that taking /i 2 = Q 2 in Eq. ( 4.1 3|) we have eliminated all dependences on 
ln(Q 2 /fi 2 ) to all orders in a. The running coupling ( 4.16 ) makes a resummation of all 
leading logarithmic corrections, i.e 



a(Q 2 



a 



n=0 



/3ia(/i 2 
2tt 



ln(Q 2 //. 2 



(4.17) 



This higher-order logarithms correspond to the contributions from an arbitrary number 
of one-loop self-energy insertions along the intermediate photon propagator in Fig. |TT 
[l-U R (q 2 /^) + (U R (q 2 /^)) 2 + •••]. 
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Figure 13: Feynman diagrams contributing to the renormalization of the strong coupling. 
The dashed loop indicates the ghost correction discussed in Appendix B. 



4.3 The QCD running coupling 

The renormalization of the QCD coupling proceeds in a similar way. Owing to the 
non-abelian character of SU(3)c, there are additional contributions involving gluon self- 
interactions. From the calculation of the relevant one-loop diagrams, shown in Fig. [T3|, 
one gets the value of the first /3-function coefficient p 



A = ^T F N f - ^C A = H^LZ_IH£ . (4 , 8) 

The positive contribution proportional to Nf is generated by the q-q loops and corresponds 
to the QED result (except for the Tp factor). The gluonic self- interactions introduce the 
additional negative contribution proportional to Nq- This second term is responsible for 
the completely different behaviour of QCD: (3\ < if Nf < 16. The corresponding QCD 
running coupling, 

u s {Q 2 ) = a ^ (g ° } > ( 4 -!9) 



2tt 



decreases at short distances, i.e. 



lim a s (Q 2 ) = 0. (4.20) 



Thus, for Nf < 16, QCD has indeed the required property of asymptotic freedom. The 
gauge self-interactions of the gluons spread out the QCD charge, generating an antiscreen- 
ing effect. This could not happen in QED, because photons do not carry electric charge. 
Only non-abelian gauge theories, where the intermediate gauge bosons are self-interacting 



particles, have this antiscreening property [24 



Although quantum effects have introduced a dependence with the energy, we still 
need a reference scale to decide when a given Q 2 can be considered large or small. An 
obvious possibility is to choose the scale at which a s enters into a strong-coupling regime 
(i.e. a s ~ 1), where perturbation theory is no longer valid. A more precise definition can 
be obtained from the solution of the /5-function differential equation Ql.14 ). At one loop, 
one gets 

ln/x+ , 2 , = InA, (4.21) 
where InA is just an integration constant. Thus, 

a - ( " 2) = -fthWAy (4 ' 22) 
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In this way, we have traded the dimensionless parameter g s by the dimensionful scale A. 
The number of QCD free parameters is the same (1 for massless quarks), but quantum 
effects have generated an energy scale. Although, Eq. ( |4.19 ) gives the impression that the 
scale-dependence of a s (/i 2 ) involves two parameters, fig and a s (/Xg), only the combination 
( 4.21 ) is actually relevant, as explicitly shown in ( 4.22|) . 

When n » A, a s (/i 2 ) — > 0, so that we recover asymptotic freedom. At lower 
energies the running coupling gets bigger; for \i — > A, a s ([i 2 ) — > oo and perturbation 
theory breaks down. The scale A indicates when the strong coupling blows up. Eq. (|4.22|) 
suggests that confinement at low energies is quite plausible in QCD; however, it does not 
provide a proof because perturbation theory is no longer valid when \i — > A. 



4.4 Higher orders 

Higher orders in perturbation theory are much more important in QCD than in 
QED, because the coupling is much bigger (at ordinary energies). Unfortunately, the 
calculations are also technically more involved. Nevertheless, many quantities have been 
already computed at 0(a 2 s ) or even 0(a^). The /3-function is known to three loops; in the 



MS scheme, the computed higher-order coefficients take the values [^5 

51 19 „ 1 

64 



02 



19 Ar 

1 N f . 

4 12 1 



03 



„„„„ 5033 Ar 
-2857 + ^—N f 



9 



325 nr2 

27 f 



(4.23) 



If Nf < 8, 02 < (03 < for Nf < 5) which further reinforces the asymptotic freedom 
behaviour. 

The scale dependence of the running coupling at higher-orders is given by: 



a s (/i 2 ) = a a (nl) 



77 



7T 



In (/i 2 //io 



+ 



(4.24) 



or, in terms of A, 



2tt 



02 



In Gu 2 /A2) 



1 (-0 1 )\n(^/A 2 ] 



In 



±ln(^ 2 /A 2 



+ 



. (4.25) 



When comparing different QCD fits to the data, it is worth while to have in mind 
that any given value of a s refers to a particular selection of scale and renormalization 
scheme. Moreover, the resulting numerical values can be different if one works at lead- 
ing (LO), next-to- leading (NLO) or next-to-next-to- leading (NNLO) order. Although the 
parameter A does not depend on the scale, it is a scheme-dependent quantity. For instance: 



A 2 



p 7b 

— A^ 

471- MS 



(4.26) 



Moreover, Alo 7^ Anlo 7^ Annlo- In fact, slightly different definitions of A can be given 
at NLO, depending on the way the integration constant is chosen when solving the 0- 
function differential equation. Moreover, since in the MS and MS schemes the /3-function 
coefficients depend on Nf, A takes different values when the number of flavours is changed. 
At NLO, the relation between the A scales for 3 and 4 flavours is given by: 



A, 



Ah 



A s 



m, 



2/25 



In (m 2 /A 



-107/1875 



(4.27) 
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Hadrons 



Figure 14: e + e — > 7*, Z* — > hadrons. 



5. PERTURB ATIVE QCD PHENOMENOLOGY 
5.1 e+e — > hadrons 

The inclusive production of hadrons in e + e~ annihilation is a good process for 
testing perturbative QCD predictions. The hadronic production occurs through the basic 
mechanism e + e~ — > 7*,Z* — ■> gg, where the final q-q pair interacts through the QCD 
forces; thus, the quarks exchange and emit gluons (and q'-q' pairs) in all possible ways. 

At high energies, where a s is small, we can use perturbative techniques to predict 
the different subprocesses: e + e~ — > qq, qqG, qqGG, . . . However, we still do not have a 
good understanding of the way quarks and gluons hadronize. Qualitatively, quarks and 
gluons are created by the q-q current at very short distances, x ~ 1/ y/s. Afterwards, 
they continue radiating additional soft gluons with smaller energies. At larger distances, 
x ~ 1/A, the interaction becomes very strong and the hadronization process occurs. 
Since we are lacking a rigorous description of the confinement mechanism, we are unable 
to provide precise predictions of the different exclusive processes, such as e + e~ — > lGn. 
However, we can make a quite accurate prediction for the total inclusive production of 
hadrons: 

a(e + e~ — > hadrons) = cr(e + e~ — > qq + qqG + qqGG + ...). (5.1) 

The details of the final hadronization are irrelevant for the inclusive sum, because the 
probability to hadronize is just one owing to our confinement assumption. 



q 



(7 ~ 




Figure 15: Diagrammatic relation between the total hadronic-production cross-section and 
the two-point function n Miy (g). The qq blob contains all possible QCD corrections. The 
dashed vertical line indicates that the blob is cut in all possible ways, so that the left 
and right sides correspond to the production amplitude T and its complex-conjugate T\ 
respectively, for a given intermediate state. 
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Well below the Z peak, the hadronic production is dominated by the 7-exchange 
contribution. Thus, we can compute the cross-sections of all subprocesses e + e~ — > 7* — > 
qq, qqG, ... (at a given order in a s ), and make the sum. Technically, it is much easier to 
compute the QCD T-product of two electromagnetic currents [J^ m = J2f QfQfl^Qf)' 



TT(q) = ijd A x (0\T (j^(x) J^(0)+) |0> = (-g^q 2 + gV) II em (g 2 



(5.2) 



As shown in Fig. IIR the absorptive part of this object (i.e. the imaginary part, which 
results from cutting -putting on shell- the propagators of the intermediate exchanged 
quarks and gluons in all possible ways) just corresponds to the sum of the squared moduli 
of the different production amplitudes. The exact relation with the total cross-section is: 



Re+e- = 



a e T e" 



hadrons) 



12tt ImlL 



(5.3) 



Neglecting the small (away from thresholds) corrections generated by the non-zero 
quark masses, the ratio R e + e - is given by a perturbative series in powers of a s (s): 



Y,Q 2 f) Ncli + Y,Fn( 

K f=l ) \ n>l V 



OCAS) 



' N 



71 



'a s (/i 2 



7T 



+ 



S 



F 3 + FAln(-]+F 1 (|ln(^)+^ln 2 (- 



s 



4 



s 



a s (jjL 



2\\ 3 



7T 



(5.4) 



The second expression, shows explicitly how the running coupling a s (s) sums an infinite 
number of higher-order logarithmic terms. 

So far, the calculation has been performed to order ac^, with the result (in the MS 



scheme) [f2~6|,|27 



Fx 
F 2 



1, 

1.986 - 0.115JV 



-6.637 - I.2OOA7 - 0.005iV 2 - 1.240 



(5.5) 



/ 



g/g/ 



Note the different charge-dependence on the last term, which is due to the contribution 
from three intermediate gluons (with a separate quark trace attached to each electromag- 
netic current in Fig. |i~5|). 

For 5 flavours, one has: 



R P 



11 I ocAs) I ocAs) 

1 + + 1.411 

3 7T V 7T 



12.80 



a.As 



71 



+ 0(a 



(5.6) 



The perturbative uncertainty of this prediction is of order af, since the coeffi- 
cient F 4 is unknown. This uncertainty also includes the ambiguities related to the choice 
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of renormalization scale and scheme. Although, the total sum of the perturbative se- 
ries is of course independent of our renormalization conventions, different choices of 
scale and/or scheme lead to slightly different numerical predictions for the truncated 
series. For instance, the perturbative series truncated at a finite order N, ~~ 



(j2fQ 2 f) Nc |l + J2n=i Fn (~^) }' nas an explicit scale dependence of order 



a N+1 : 



dR 



W /„ t,.2\\ N+1 



e T e" 




The numerical values of a s and the F n (n > 2) coefficients depend on our choice of 
scheme (also (3 n for n > 3). For instance, at second ordeiQ, the relation between the MS 
and MS schemes is: 



8 n> MS , 

l + ^[ln(47r)- 7£ ]^ + ...^ (5.8) 



f ms = f ms _ F & [ln ^ _ = 7>35g _ _4 41iV/ _ (5 9) 

The difference between both schemes is obviously a higher-order effect. With Nf = 5, the 
MS scheme leads to a second-order coefficient F 2 MS = 5.156, which is a factor 3.6 bigger 
than F 2 MS . Thus, the perturbative series looks more convergent with the MS choice. 

The theoretical prediction for R e + e - (s) above the b-b threshold is compared [JTOj] in 
Fig. § with the measured data, taking into account mass-corrections and electroweak {Z- 
exchange) contributions. The two curves correspond to A^ _5 ' ) = 60 MeV (lower curve) 
and 250 MeV (upper curve). The rising at large energies is due to the tail of the Z peak. 
A global fit to all data between 20 and 65 GeV yields [p8]| 

« s (34GeV) = 0.146 ±0.030. (5.10) 

The hadronic width of the Z boson can be analyzed in the same way: 



The global factor 



)EW _ T,f(v} + a 2 f ) 



Rz = ^T—^ (1 + *ew) (5-12) 

contains the underlying electroweak Z — > Y^fQfQf decay amplitude. Since both vector and 
axial- vector couplings are present, the QCD-correction coefficients F n are slightly different 
from F n for n > 2. For instance, the Z axial coupling generates the two-loop contribution 
Z — » tt —>■ GG — > qq (through triangular quark diagrams), which is absent in the vector 
case; this leads to an additional O (a 2 mf/M^) correction. 

In order to determine a s from Rz, one needs to perform a global analysis of the 
LEP/SLC data, taking properly into account the higher-order electroweak corrections 
|]2"9l, [3011 . The latest a s value reported by the LEP Electroweak Working Group [^TJ is 

a s (M 2 z ) = 0.125 ±0.005 ±0.002. (5.13) 

2) Actually, at second order a scheme is completely specified by a single parameter. Thus, scale 
and scheme dependence is just the same at this order. The relations in Eqs. ( fTsj ) and (5.9) 



are equivalent to a change of scale: /U^g = (47r/e 7E ) \ih 



MS 
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5.2 t — ► u T + hadrons 

The calculation of QCD corrections to the inclusive decay of the r lepton P2"|-p(3 
looks quite similar from a diagrammatic point of view. One just puts all possible gluon 
(and qq) corrections to the basic decay diagram in Fig. [|, and computes the sum 

T(t~ — ► v T + hadrons) = r(r~ — ► z/ T + qq) + T(t~ — > z/ T + ggG) + r(r~ — > z/ T + qqGG) + ■ ■ ■ 

(5.14) 

As in the e + e~ case, the calculation is more efficiently performed through the two-point- 
function 

nf(g) = • / ^e" 1 (0|T (^(x)L^(O)t) |0) = (-g^q 2 + <f<f) Up \q 2 ) + q^U?\q 2 ) , 

(5.15) 

which involves the T-ordered product of two left-handed currents, L M = wy^{l — rf)do. 
This object can be easily visualized through a diagram analogous to Fig. [T3|, where the 
photon is replaced by a W~ line and one has a tv t pair in the external fermionic lines 
instead of the e + e~ pair. The precise relation with the ratio R T is: 

r(r~ — > z/ r + hadrons) 



I — > z/ r e z/ e ) 



12tt 




1 + 2^7 ) Imni 1} (s) + Imni 0) ' 

mi 



(5.16) 



The three-body character of the basic decay mechanism, r~ — > v T udg, shows here a 
crucial difference with e + e~ annihilation. One needs to integrate over all possible neutrino 
energies or, equivalently, over all possible values of the total hadronic invariant-mass s. 
The spectral functions ImlT^^s) contain the dynamical information on the invariant- 
mass distribution of the final hadrons. The lower integration limit corresponds to the 
threshold for hadronic production, i.e. m n (equal to zero for massless quarks). Clearly, 
this lies deep into the non-perturbative region where confinement is crucial. Thus, it is 
very difficult to make a reliable prediction for the integrand in Q5.16 ). 



Im(s) 




Re(s) 



Figure 16: Integration contour in the complex s-plane used to obtain Eq. Q5.17Q . 

Fortunately, we have precious exact (i.e. non-perturbative) information on the 
dynamical functions lTj° (s), which allows us to accurately predict the total integral 
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( |5.16| ): n^ 0,1 ^(s) are analytic functions in the complex s-plane except for a cut in the 
positive real axis. The physics we are interested in lies of course in the singular region, 
where hadrons are produced. We need to know the integral along the physical cut of 



mill 



(0,1), 



s) 



(0,1), 



s + it) 



it)}. However, we can use Cauchy's theorem 



(close integrals of analytic functions are zero if there are no singularities within the in- 
tegration contour), to express R T as a contour integral in the complex s-plane running 
counter-clockwise around the circle \s\ — m 2 fB2T|34 



i? x = Qui 



J\s\=m% m$ 



mi 




mi 



(5.17) 



The advantage of this expression is that it requires dynamical information only for complex 
s of order m 2 , which is significantly larger than the scale associated with non-perturbative 
effects in QCD. A perturbative calculation of R T is then possible. 

Using the so-called Operator Product Expansion techniques it is possible to show 
33]. |34], |36| that non-perturbative contributions are very suppressed [~ (A/m r ) 6 ]. Thus, 



R T is a perfect observable for determining the strong coupling. In fact, r decay is prob- 
ably the lowest energy process from which the running coupling constant can be ex- 
tracted cleanly, without hopeless complications from non-perturbative effects. The r mass, 

|10| , lies fortuitously in a compromise region where the coupling 
is large enough that R T is very sensitive to its value, yet still small enough 
that the perturbative expansion still converges well. 



m 7 
constant a 



1.7771 ±8:Soot GeV 



The explicit calculation gives |34|, |36[1 : 



Rr = 3 (\V ud \ 2 + \V US \ 2 ) Sew {l + S' EW + 5<°> + <5 NP } 



(5.1* 



where Sew = 1-0194 and S' EW = 0.0010 are the leading and next-to-leading electroweak 
corrections, and 5^ contains the dominant perturbative-QCD contribution: 



5 (o) 



a Am: 



+ 



7T 
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+ 

a s (m 2 T ) 



FL - — F 2 fr - — fa 
3 12 24 



+ 5.2023 



'a s {m 2 T ) 



71 



71 




+ 0(at 



+ 0(a 4 s ) . 



(5.19) 
(5.20) 



The remaining factor <5np ~ —0.016 ± 0.005 includes the estimated |34],|3()| small mass- 
corrections and non-perturbative contributions. 

Owing to its high sensitivity to a s [[33|,[34[ the ratio R r has been a subject of intensive 
study in recent years. Many different sources of possible perturbative and non-perturbative 
contributions have been analyzed in detail. Higher-order logarithmic corrections have been 
resummed ]35| , leading to very small renormaliz at ion- scheme dependences. The size of the 
non-perturbative contributions has been experimentally analyzed, through a study of the 
invariant-mass distribution of the final hadrons ||37|| ; the present data implies J38| 5 



NP 



(0.3 ±0.5)% confirming the predicted [34]] suppression of non-perturbative corrections. An 
exhaustive summary of the R T analysis can be found in Ref. [pH]. 
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Using the Particle Data Group values for the r lifetime and leptonic branching 
ratios [TTJ, the theoretical analysis of R T results in a fitted value of a s 

(5.21) 



0.33 ±0.03 



which is significantly larger than (|5.13|) . After evolution up to the scale Mz, the strong 



coupling constant in Eq. (|5.21|) decreases to a s (M|) = 0.120io]ooi) m excellent agreement 
with the Z-width determination and with a smaller error bar. This comparison provides 
a beautiful test of the predicted running of a s . 



5.3 



e + e 



jets 





Figure 17: Gluon bremsstrahlung corrections to e + e 



qq. 



At lowest-order in the strong coupling, the hadronic production in e + e~ collisions 
proceeds through e + e~ — > qq. Thus, at high-energies, the final hadronic states are pre- 
dicted to have mainly a two-jet structure, which agrees with the empirical observations. 
At 0(a s ), the emission of a hard gluon from a quark leg generates the e + e _ — » qqG tran- 
sition, leading to 3-jet configurations. For massless quarks, the differential distribution of 
the 3-body final state is given by: 



1 d 2 a 

<t dx\dx- A 



2a. 



rf ^ I ,y ^~ 



37T (1 — Xl)(l — X 2 ) 



where 
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Aira 2 
~3s~ 



(5.22) 



(5.23) 



qq cross-section. The kinematics is defined through the 

^2 



is the lowest-order e + e~ — > 7* 

invariants s = q 2 and Sjj = (pj + p^) 2 = (g — pk) 2 = s(l — x^) (z, j, = 1, 2, 3), where pi, 
p 2 and p 3 are the quark, antiquark and gluon momenta, respectively, and q is the total 
e + e~ momentum. For given s, there are only two independent kinematical variables since 

x 1 + x 2 + x 3 = 2. (5.24) 

In the centre-of-mass system [g M = (y/s,Q)], x, = Ei/E e = lEij^fs. 

Eq. ([5.22 ) diverges as ii or x 2 tend to 1. This is a very different infinity from 
the ultraviolet ones encountered before in the loop integrals. In the present case, the tree 
amplitude itself is becoming singular in the phase-space boundary. The problem originates 
in the infrared behaviour of the intermediate quark propagators: 

Xt -»• 1 (p 2 +P3) 2 = 2 (p 2 -Pz) -> ; 

x 2 -> 1 (pi + p 3 ) 2 = 2 (p x • ps) -»• . 

There are two distinct kinematical configurations leading to infrared divergences: 



(5.25) 
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1. 



Collinear gluon: The 4-momentum of the gluon is parallel to that of either the 
quark or the antiquark. This is also called a mass singularity, since the divergence 
would be absent if either the gluon or the quark had a mass (p3 ||p2 implies S23 = if 



vl = V% = 0). 



2. 



Soft gluon: p 3 — > 0. 

In either case, the observed final hadrons will be detected as a 2-jet configuration, 
because the qG or qG system cannot be resolved. Owing to the finite resolution of any 
detector, it is not possible (not even in principle) to separate those 2-jet events generated 
by the basic e + e~ — > qq process, from e + e~ — > qqG events with a collinear or soft gluon. 
In order to resolve a 3-jet event, the gluon should have an energy and opening angle (with 
respect to the quark or antiquark) bigger than the detector resolution. The observable 
3-jet cross-section will never include the problematic region x\,2 —* 1; thus, it will be finite, 
although its value will depend on the detector resolution and/or the precise definition of 
jet (i.e. a depends on the chosen integration limits). 

On the other side, the 2-jet configurations will include both e + e~ — > qq and 
e + e~ — > qqG with an unobserved gluon. The important question is then the infrared 
behaviour of the sum of both amplitudes. The exchange of virtual gluons among the 
quarks generate an 0(a s ) correction to the e + e~ — > qq amplitude: 



T[e H 



qq] =T + T 1 + 



(5.26) 



where T is the lowest-order (tree- level) contribution, T x the 0(a s ) correction, and so 
on. The interference of T and T\ gives rise to an 0(a s ) contribution to the e + e~ — > qq 
cross-section. 




+ 




+ 




Figure 18: 1-loop gluonic corrections to e + e 



qq. 



We know already that loop diagrams have ultraviolet divergences which must be 
renormalized. In addition, they also have infrared divergences associated with collinear 
and soft configurations of the virtual gluon. One can explicitly check that the 0(a s ) 
infrared divergence of a(e + e~ — > qq) exactly cancels the one in a(e + e~ — > qqG), so that 
the sum is well-defined: 



cr(e + e — > qq) + cr(e + e — ► qqG) 



1 + — + 

7T 



(5.27) 



This is precisely the inclusive result discussed in Sect. |5.1| . This remarkable cancellation of 
infrared divergences is actually a general result (Bloch-Nordsieck |M| and Kinoshita-Lee- 
Nauenberg f|D| theorems): for inclusive enough cross-sections both the soft and collinear 
infrared divergences cancel. 

While the total hadronic cross-section is unambiguously defined, we need a precise 
definition of jet in order to classify a given event as a 2-, 3-, . . . , or n-jet configuration. 
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Figure 19: 2-jet configuration. 



Such a definition should be free of infrared singularities, and insensitive to the details of 
the non-perturbative fragmentation into hadrons. A popular example of jet definition is 



the so-called JADE algorithm which makes use of an invariant-mass cut y: 



3 jet <{=> Si , = ( Pl + Pj ) 2 >ys = 1,2,3) . (5.28) 

Clearly, both the theoretical predictions and the experimental measurements de- 
pend on the adopted jet definition. With the JADE algorithm, the fraction of 3-jet events 
is predicted to be: 



where 

Li a (z) = - f -^-ln£. (5.30) 
Jo 1-4 

The corresponding fraction of 2-jet events is given by R<± = 1 — -R3. The fraction of 2- or 
3-jet events obviously depends on the chosen cut y. The infrared singularities are manifest 
in the divergent behaviour of R3 for y — > 0. 

At higher-orders in a s one needs to define the different multi-jet fractions. For 
instance, one can easily generalize the JADE algorithm an classify a {pi,p%, ■ ■ ■ ,p n } event 
as a n-jet configuration provided that s^- > ys for all i,j = 1, . . . , n. If a pair of momenta 
does not satisfy this constraint, they are combined into a single momentum and the event 
is considered as a (n — 1) jet configuration (if the constraint is satisfied by all other 
combinations of momenta). The general expression for the fraction of n-jet events takes 
the form: 

K(s.y) = (£^)"" , £cf>( v ) (^flj , (5.31) 

With En-Rn = 1- 

A few remarks are in order here: 

— The jet fractions have a high sensitivity to a s [R n ~ «™~ 2 ]- Although the sensitivity 
increases with n, the number of events decreases with the jet multiplicity. 

— Higher-order a s {^ 2 y In (s//x 2 ) terms have been summed into a s {s). However, the co- 
efficients Cj(y) still contain ln k (y) terms. At low values of y, the infrared divergence 
(y — > 0) reappears and the perturbative series becomes unreliable. For large y, the 
jet fractions R n with n > 3 are small. 
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— Experiments measure hadrons rather than partons. Therefore, since these observables 
are not fully inclusive, there is an unavoidable dependence on the non-perturbative 
fragmentation into hadrons. This is usually modelled through Monte Carlo analyses, 
and introduces theoretical uncertainties which need to be estimated. 
Many different jet algorithms and jet variables (jet rates, event shapes, energy 

correlations, . . . ) have been introduced to optimize the perturbative analysis. In some 

cases, a resummation of a s (s) n ln m (y) contributions with m > n has been performed to 

improve the predictions at low y values |42|| . 
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Figure 20: Energy dependence of 3-jet event production rates Rz{y = 0.8), compared 
with predictions of analytic 0{a 2 s ) QCD calculations, with the hypothesis of an energy 
independent a s and with the abelian vector theory in 0(a\) (taken from Ref. [|43|j ). 



Fig. |20] |43| shows the energy dependence of the measured 3-jet production fraction 
i?3 (y = 0.08), compared with QCD predictions. The data is in good agreement with QCD 
and fits very well the predicted energy-dependence of the running coupling. A constant 
value of a s cannot describe the observed production rates. The figure shows also the 
predictions obtained with an abelian vector theory at 0(a\), which are clearly excluded. 

Several measurements of a s , using different jet variables, have been performed. 
All measurements are in good agreement, providing a good consistency test of the QCD 
predictions. Combining the results from all experiments at LEP and SLC, one gets the 
average value |7[] : 



0.119 ±0.006 
0.123 ±0.006 



(resummed calculations) 



(5.32) 



The two numbers correspond to different theoretical approximations used in the fits to 
extract a s . 

3-jet events can also be used to test the gluon spin. For a spin-0 gluon, the dif- 
ferential distribution is still given by Eq. ( |5.22|) , but changing the x\ + x\ factor in the 
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numerator to x|/4. In general, one cannot readily be sure which hadronic jet emerges via 
fragmentation from a quark (or antiquark), and which from a gluon. Therefore, one adopts 
instead a jet ordering, X\ > xi > 23, where x\ refers to the most energetic jet, 2 to the 
next and 3 to the least energetic one, which most likely would correspond to the gluon. 
When X2 — > 1 (xi —* 1) the vector-gluon distribution is singular, but the corresponding 
scalar-gluon distribution is not because at that point £3 = (1 — x\) + (1 — £2) — > 0. The 
measured distribution agrees very well with the QCD predictions with a spin-1 gluon; a 
scalar gluon is clearly excluded. 
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Figure 21: 68% and 95% CL contours 
in the Tf/Cf versus Ca/Cf plane, from 
OPAL data |[44|| . Expectations from vari- 
ous gauge models are also shown. 
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Figure 22: Summary of colour-factor mea- 
surements ||. The results refer to 5 active 
flavours with Tr = NfTp = 5Tp. 



The predictions for jet distributions and event shapes are functions of the colour- 
group factors T F = 1/2, Cp = (Nq — 1)/(2Nc) and Ca = N c . These quantities, defined in 
Eq. ( |A.5| ), result from the colour algebra associated with the different interaction vertices, 
and characterize the colour-symmetry group. If the strong interactions were based on a 
different gauge group, the resulting predictions would differ in the values of these three 
factors. Since the vertices contribute in a different way to different observables, these 



colour factors can be measured by performing a combined fit to the data. Fig. |21] compares 
a recent OPAL determination jyj of Ca/Cf and Tf/Cf with the values of these two ratios 
for different colour groups. The data is in excellent agreement with the SU (3) values, and 
rules out the Abelian model and many classical Lie groups. Notice that those groups shown 
by the open squares and circles are already excluded because they do not contain three 
colour degrees of freedom for quarks. Similar results have been presented by the other LEP 
experiments (the hadronic production of jets at pp colliders has also been analyzed in a 
similar way). A summary of the colour-factor ratios obtained by the different experiments 



is given in Fig. ^2. 
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6. DEEP INELASTIC SCATTERING 



We saw in Section [L^ how the deep inelastic scattering (DIS) e"p — > e~X can be 
used to learn about the proton structure. Since this involves a bound hadronic state -the 
proton-, non-perturbative phenomena such as confinement plays here a crucial role. At 
the same time, the data obeys Bjorken scaling which manifests the asymptotic freedom 
property of the strong interactions. Thus, DIS appears to be an interesting place where 
to investigate both perturbative and non-perturbative aspects of QCD. 

DIS can be visualized as a two-step process. First, the hard intermediate photon, 
which is far off its mass-shell, scatters off a quark or gluon with a large momentum transfer; 
this scattering can be adequately described by perturbation theory. Second, the outgoing 
partons recombine into hadrons in a time of 0(1/ A). Although this recombination is not 
calculable in perturbation theory, the details of the non-perturbative hadronization can 
be avoided, by considering fully inclusive rates, so that perturbative QCD can be applied. 
However, the hadronic bound-structure of the initial proton state, still introduces a non- 
perturbative ingredient: the proton structure functions. 

6.1 Free parton model 

Let us ignore any QCD interactions and let us assume that the nucleon (either 
proton or neutron) constituents are free spin-| partons. Within the quark model, the 
nucleons have three point-like constituents (p = u v u v d v , n = u v d v d v ), which we will call 
valence quarks. Gluons are of course there; however, they do not interact directly with 
the photon probe. The photon-gluon interaction only occurs through the virtual q-q pairs 
coupled to the gluon constituents. Thus, instead of gluons, the photon feels a sea of q-q 
partons within the nucleon. 

Let us denote u(x), u(x), d(x), d(x), s(x), s(x), . . . the probability distributions 
for u, u, d, d, s, s, . . . quarks with momentum fraction x in the proton. We have seen in 



Section [L2| that, within the parton model, the proton structure functions have a simple 



form in terms of parton distributions: 

F* p (x)/x = 2Fl P (x) = ^ \u(x) + u(x)\ + \ \d(x) + d(x)\ + \ \s(x) + s(x)\ + ■■■ (6.1) 

9 9 9 

The same parton distributions occur in other DIS processes such as up — > l~X or up — > 
l + X. However, since the quark couplings of the intermediate bosonic probe (a W ± in that 
case) are not the same, different combinations of these functions are measured: 

F% p (x)/x = 2F? p (x) = 2 [d(x) + s(x) + u(x) + c(x) + ■••], 
F^ p (x) = 2 [d(x) + s(x) - u(x) - c(x) + ■■■], 

(6.2) 

Fl v (£)/x = 2Fl P (x) = 2 \u(x) + c(x) + d(x) + s(x) + •••], 
F* p (x) = 2 [u(x) + c(x) - d(x) - s(x) + •■■]. 

Using isospin symmetry, we can further relate the up- and down-quark distributions in a 
neutron to the ones in a proton: 

u n ( x ) = d p (x) = d(x) ; d n (x) = u p (x) = u(x) ; (6.3) 
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the remaining parton distributions being obviously the same. Thus, combining data from 
different DIS processes, it is possible to obtain separate information on the individual 
parton distribution functions. 

The quark distributions must satisfy some constraints. Since both the proton and 
the neutron have zero strangeness, 

1 

dx [s(x) - s(x)] = . (6.4) 

o 

Similar relations follow for the heavier flavours (c, . . . ). The proton and neutron electric 
charges imply two additional sum rules, 

dx [u(x) — u(x)] = 2 , dx [d(x) — d(x)\ = 1 , (6.5) 

o Jo 

which just give the excess of u and d quarks over antiquarks. 

The quark-model concept of valence quarks gives further insight into the nucleon 
structure. We can decompose the u and d distribution functions into the sum of valence 
and sea contributions, and take the remaining parton distributions to be pure sea. Since 
gluons are flavour singlet, one expects the sea to be flavour independent. In this way, the 
number of independent distributions is reduced to three: 

u(x) = u v (x) + q s (x) , 

d(x) = d v (x) + q s (x) , (6.6) 
u(x) = d(x) = s(x) = s(x) = . . . = q s (x) . 



Within this model, the strangeness sum rule ( |6.4| ) is automatically satisfied, while 
imply constraints on the valence-quark distributions alone. 

In the analogous situation of quasi-elastic electron-deuterium scattering, the ob- 
served structure function shows a narrow peak around x = \. This is to be expected, 
since the deuteron has two nucleon constituents with « \Md which share the total 
momentum in equal terms. A simple three-quark model for the nucleon would suggest 
the existence of a similar peak at x — | in the proton and neutron structure functions. 
However, the distribution shown in Fig. [7| does not show such behaviour. The difference 
can be easily understood as originating from the parton-sea contributions. Taking the 
difference between the proton and neutron structure functions, where the contribution 
from the sea cancels, the data exhibits indeed a broad peak around x = |. 



Our isospin symmetric parton model implies the so-called Gottfried sum rule ||45|| : 

C — [F?(x) - F?(x)] = \fdx K(x) - d v {x)\ = ~ , (6.7) 
jo x 6 Jo o 

which is well satisfied by the data. Another interesting quantity is the ratio 

F|" (x) 4:d v (x) + uJx) + E sea 



F 2 ep (x) 4u v (x) + d v (x) + S 5 



(6. 



where S sea is the total sea contribution. Since all probability distributions must be 
positive-definite, this ratio should satisfy the bounds \ < F 2 e ™ \x) / F 2 ep \x) < 4, which are 
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consistent with the data. The measured ratio appears to tend to 1 at small x, indicating 
that the sea contributions dominate in that region. 

The conservation of the total proton momentum implies an important sum rule: 



u(x) + u(x) + d(x) + d(x) + s(x) + s(x) + 



1 -e 



(6.9) 



where e is the fraction of momentum that is not carried by quarks. One finds experimen- 
tally that e pa | (at Q 2 ~ 10-40 GeV 2 ), suggesting that about half of the momentum is 
carried by gluons. This shows the important role of gluons in the proton structure. Al- 
though the naive quark model works very well in many cases, it is a too gross simplification 
as a model of hadrons, at least al large Q 2 . 



6.2 QCD-improved parton model 

At lowest order the DIS process occurs through the hard scattering between the 
virtual photon (W, or Z) and one constituent parton. The obvious first QCD corrections 
will be due to real gluon emission by either the initial or final quark. To get rid of infrared 
divergences, the one- loop virtual gluon contribution should also be taken into account. 




Small Q 2 >► Large Q 2 

Figure 23: Resolution of the photon probe as function of Q 2 . 

One can easily understand the main qualitative features of gluon emission, with 
a few kinematical considerations. At very low values of momentum transfer, the proton 
behaves as a single object, either point-like (at Q 2 ~ 0), or with a finite size. At higher 
energies, the photon is sensitive to shorter distances and scatters with the constituent 
partons. Increasing further the momentum transfer, the photon probe has a greater sen- 
sitivity to smaller distances, and it is able to resolve the scattered quark into a quark and 
a gluon. Thus, a parton with momentum fraction x can be resolved into a parton and a 
gluon of smaller momentum fractions, x' < x and x — x', respectively. In a similar way, a 
gluon with momentum fraction x can be resolved into a quark and an antiquark. 

This simple picture implies that increasing the Q 2 , the photon will notice some 
qualitative changes in the parton distributions: 

— Gluon bremsstrahlung will shift the valence and sea distributions to smaller x values. 

— The splitting of a gluon into a quark-antiquark pair will increase the amount of sea 
(mostly at small x). 

Thus, without any detailed calculation, one can expect to find a definite Q 2 dependence 
in the parton distributions; i.e. violations of Bjorken scaling due to the underlying QCD 
interactions. 



31 





+ >>, . — + >> - i > + 




Figure 24: Leading gluonic correction to the basic DIS parton process. 

Let us consider a quark with momentum fraction y. At lowest order, its contribution 
to the proton structure function can be written as 

2F ( ?\x) = e\ f*dyq(y)5(x-y) . (6.10) 

If the quark emits a gluon before being struck by the photon, its momentum fraction 
will be degraded to yz (0 < z < 1). Assuming that the quark remains approximately 
on-shell, (q + yzP) 2 m m 2 m 0, implying that yz = Q 2 /2(P ■ q) = x. Therefore, Fi(x) 
gets contributions from quarks with initial momentum fractions y > x. 

The explicit calculation of the diagrams in Fig. [24] gives the result: 

2AFi q \x) = e 2 q ^- £dyq{y) £ dz5(yz-x) {P q + q (z) \n(Q 2 /u 2 m ) + C(z)} , (6.11) 



where 



P qq (z) ee C F (\±4) ( 6 - 12 ) 



is called the quark splitting function. 

The important feature in Eq. ( |6.11| ) is the appearance of a scaling violation through 
the logarithmic a s correction. A careful analysis of the different Feynman diagrams shows 
that ultraviolet divergences are absent in the total contribution. Therefore, this loga- 
rithm has a completely different origin than the ultraviolet ones found in Section ^] The 
logarithmic behaviour is now generated by infrared singularities of the type discussed 
in Section |5.3|. More precisely, there is a collinear singularity associated with the gluon 



emission process, which has been regulated with the infrared cut-off z/rR. 

The general theorems on the cancellation of infrared divergences do not protect the 
structure function Fi (x), because this quantity is not inclusive enough. The divergence 
shows up when one tries to resolve the original quark with momentum fraction y into 
a quark with momentum fraction yz and a gluon. P qq {z) is just the coefficient of the 
logarithmic divergence associated with the splitting process q — > qG. Physical observables 
should not depend on any cut-off, however, our definition of a parton distribution obviously 
depends on the power resolution of our photon probe. While at low Q 2 the photon was 
testing a single parton with momentum fraction y = x, now it feels the splitting of a 
quark with y > x into a quark and a gluon with separate parton distributions. 

The divergence should then be reabsorbed into the observable parton distribution 
function: 

q(x, Q 2 ) = q(x, u 2 m ) + %L HQ 2 /v\ K ) f *H q(y) P+ (x/y) . (6.13) 

ait j x y 
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Both the bare distribution q(x, i/f R ) and the a s correction depend on the infrared cut- 
off, but this dependence cancels out and does not show up in the physical distribution 
functionQg(x, Q 2 ). Instead, the parton distribution is now a Q 2 -dependent quantity, which 
fits with our intuitive picture that the photon probe increases its resolution power with 
the scale. In terms of q(x, Q 2 ), the contribution of the quark q to the proton structure 
function is given by: 

2F®(x) = e 2 q \q(x,Q 2 ) + ^ f *H q(y) C (x/y)\ . (6.14) 
[ 2tt Jx y J 

The individual diagrams in Fig. [24] have also a soft-gluon singularity, which mani- 
fests in the divergent behaviour of P qq (z) at z = 1. This singularity cancels exactly in the 
total sum of the gluon-emission and virtual-gluon-exchange contributions. The net result 
is a slight modification in the definition of the splitting function: 

P+(z) 5(yz - x) = P qq (z) [8(yz - x) - 5{y - x)\ . (6.15) 

Eq. ( |6.13| ) shows an important thing: although perturbative QCD is not able to 
predict the actual value of the distribution function, it does predict how this distribution 
evolves in hi(Q 2 ). Thus, given its value at some reference point Qg, one can compute 
the quark distribution at any other value of Q 2 (high-enough for perturbation theory 
to be valid). Including the leading higher-order logarithmic corrections into the running 
coupling, the (Revolution of the parton distribution is given by [[46|,|4"7||: 

q2 j& q(x ' q2) = r f q{y ' q2) p ^ x/y) • (6 - i6) 

Thus, the change in the distribution for a quark with momentum fraction x, which in- 
teracts with the virtual photon, is given by the integral over y of the corresponding 
distribution for a quark with momentum fraction y > x which, having radiated a gluon, is 
left with a fraction x/y of its original momentum. The splitting function has then a very 
intuitive physical interpretation: (a s /2ir)P+(x/y) is the probability associated with the 
splitting process q(y) — > q(x)G. This probability is high for large momentum fractions; 
i.e. high-momentum quarks lose momentum by radiating gluons. Therefore, increasing 
Q 2 , the quark distribution function will decrease at large x and will increase at small x. 

(X-X')P ( X " X ')P , , ND 

1 ' (x-x)P 

x'P xP 'P xP H^jx 'P 

Figure 25: Basic parton-splitting processes. 





3) Notice, however, that the precise definition of q(x,Q 2 ) is factorization-scheme dependent, 
since we could always include some arbitrary non-logarithmic a s correction into q(x,Q 2 ), by 
simply shifting the C (x/y) correction factor in fl6.140 . 
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The evolution equation ( |6 . 1 is only correct for non-singlet distributions such as 
qi(x) — qj(x), where the (flavour-singlet) gluon contribution cancels out. In general, one 
needs also to consider the effects coming from the splitting of a gluon into a quark and 
an antiquark, which interact with the photon probe. The obvious generalization is p7| : 



Q 2 



d 



q(x,Q 2 ) 
dQ 2 I G(x,Q 2 ) 



a s (Q 2 ) f l dy 
2it Jx y 



P+(x/y) P qG {x/y) 
P Gq (x/y) P GG (x/y) 



q(y,Q 2 ) 
G(y,Q 2 ) 



(6.17) 



where Pa q (z) = P qq {l — z) determines the probability that a quark radiates a gluon with 
a fraction z of the original quark momentum, while 



P qG {z) 
Pgg(z) 



2C A 



z z + (l-z 

z 



(6.18) 



1 — z 
+ + z{l - z) 



1 



+ ~(llC A -W f T F )8(l-z) , (6.19) 
6 



are the gluon splitting functions into qq and GG, respectively. The subindex "+" in the 
1/(1 — z) + factor indicates that the z — 1 divergence disappears through 



dzf{z)[g{z)\. 



dz[f(z)-f(l)]g(z) . 



(6.20) 



6.3 Moments of the structure functions 

The previous discussion has been based on rather qualitative arguments. Never- 
theless, the predicted evolution equation can be derived on a more rigorous basis using 
the formal framework of the operator product expansion [4~8|, which allows to make a full 
QCD analysis of the moments 

M q N {Q 2 ) = ( l dxx N - l q{x,Q 2 ) ; M%(Q 2 ) = [' dx x N ~ x G(x, Q 2 ) . (6.21) 
Jo Jo 



Taking moments on both sides of Eq. ( p. 17 ), one finds 
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where 



7^ = I dzz 



(6.22) 



(6.23) 



are the so-called anomalous dimensions. Performing the trivial integrals, one gets: 

1 



Iqq 



Iql 



l N Gq 



iGG = 2C A 



c F 

C F 



1 

2 



N{N +1) 
2 + N + N 2 
N(N + l)(JV + 2) 
2 + + A^ 2 
N{N 2 - 1) ' 



N i 
k=2 K . 



(6.24) 



1 



+ 



1 



iV 



+ 



12 N(N-l) (N+l)(N + 2) 



-E 

k=2 
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For a non-singlet structure function, where the gluon component is absent, the 
evolution differential equation leads to the solution 

Mr m = M r( Qi) ggg) " ; d „ s 7 » /ft = ^ . 

The first moment has d\ — 0; therefore, the Gottfried sum rule ( |6.7| ) does not get any 
QCD correction at this leading order. For iV > 2, djy > so that M^ ns (Q 2 ) decreases as 
Q 2 increases, indicating a degradation of momentum in the non-singlet quark distribution. 

Let us now consider the flavour-singlet structure function = 5Di [<li( x ) + Qi( x )}- 
The N = 2 moments Mf(Q 2 ) and M§(Q 2 ) give the average total fraction of momentum 
carried by quarks and gluons, respectively. The corresponding coupled evolution equations 
can be easily solved. The sum of both moments does not depend on Q 2 , since the total 
momentum is conserved: 

M?(Q 2 ) + M?(Q 2 ) = 1. (6.26) 
The evolution of the N = 2 singlet distribution then takes the simple form 



Mf{Q 2 ) - §- F M?{Q 2 ) M 2 S (Q 2 ) - TeS^ (* s {Q 2 )\ d > 
MKQl)-&M§{Ql) Mf(Q2)__^_ WqD) 



(6.27) 

Ic^"2 v^o; wo; - ie+3^ 1 n " 7 "' ' 

with cf = 2{AC F + iV/)/(33 - 2iV/). If iV> < 16, cf > and the right-hand side will 
decrease for increasing Q 2 . Thus, one gets a prediction for the asymptotic values of the 
average total momentum carried by quarks and gluons: 

lim M 2 S (Q 2 ) = ; Hm M 2 G (Q 2 ) = -i^— . (6.28) 

Q 2 ^oo 16 + SJMf Q 2 ^oo 16 + 3A/j 

For Nf = 4, this gives | and |, in good agreement with the empirical observation that 
for Q 2 in the range 10-40 GeV 2 each fraction is very close to |. 

A very interesting issue is the behaviour of the parton distributions at the end- 
points x = and x = 1. The large iV moments probe the x — > 1 region, while the 
low x behaviour is controlled by the iV — ► 1 limit. As iV increases, 7^ and 7q 9 tend 
to zero, so that the evolution equations (|6.22|) decouple; i.e. the large x behaviour of the 



quarks is independent of the gluon evolution. When x — > 1 the gluon distribution function 
approach zero more rapidly than the quark ones. For large values of x the quark content 
of the nucleon is the relevant one. Notice that x — 1 means W 2 = M 2 , i.e. it actually 
corresponds to the elastic photon-nucleon scattering. 

At low x, x/y — > and the splitting functions Pgg{ x /u) an d PG q {x/y) diverge. The 
gluon distribution function becomes then dominant. The low x behaviour is controlled by 
the singular iV — > 1 limit of the gluon anomalous dimension 7g G ~ 2Ca/{N — 1). Making 
a saddle-point approximation, the iV — > 1 moment can be inverted; one finds in this way 
that for low x the gluon distribution function behaves as 




G{x) expWC(Q 2 ) In - , (6.29) 



with C(Q 2 ) a calculable function. Obviously, this behaviour cannot be true for arbitrarily 
small x; something must stop the growing of the gluon distribution before running into 
unitarity problems. 
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Figure 26: Recent measurements of F 2 ep (x, Q 2 ) p9|-|52|. The F 2 ep (x, Q 2 ) values are plotted 
with all but normalization errors in a linear scale adding a term c(x) = 0.6 (i x — 0.4) to 
F2, where i x is the bin number starting at i x — 1 for x = 0.13. The curves represent a 
phenomenological fit to the data. (Taken from Ref. f49|j). 



Kinematically, low x means the high-energy (high W 2 ) limit for the virtual photon- 
nucleon scattering. The e-p HERA collider is ideally suited for studying this region. The 
HERA experiments extend the previously accessible kinematic range up to very large 
squared momentum transfers, Q 2 > 5 x 10 4 GeV 2 , and to very small values of x < 10 -4 . 



The measurements reported so far [49],[5(| observe indeed a significant rise of the structure 
function -F 2 ep (x, Q 2 ) with decreasing x, at fixed Q 2 . Around x ~ 10~ 3 the decrease of x 
by an order of magnitude amounts to a rise of F^ix, Q 2 ) of about a factor of two. The 
observed Q 2 behaviour is consistent with the expected scaling violations, i.e. a weak rise 
of F2 P (x, Q 2 ) with increasing Q 2 for x < 0.1. The most recent data [49-52 on the proton 
structure function F^lx, Q 2 ) are shown in Figs. E6| and |27. 



6.4 QCD fits to DIS data 

There is a twofold motivation for making careful analyses of DIS data. First, the 
experimental measurement of the parton distributions provides very valuable informa- 
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Figure 27: x dependence of the measured structure function F| p (x, Q 2 ), for different Q 2 
values [19, ET), p>2]] . The curves represent a phenomenological fit to the data. (Taken from 
Ref. M- 



tion on the non-perturbative regime of the strong interactions (in addition, these parton 
distributions are needed for making predictions of hard-scattering processes in hadronic 
collisions). Second, the measured Q 2 evolution (the slopes of the distributions) can be 
compared with perturbative QCD predictions. 

Usually, one adopts some motivated parametrization of the quark and gluon distri- 
butions at a fixed momentum-transfer Qq. The evolution equations are then used to get 
the proton (or neutron) structure functions at arbitrary values of Q 2 , and a global fit to 
the data is performed. 

In the actual analysis one needs to worry about the unavoidable presence of ad- 
ditional non-perturbative contributions. The perturbative evolution equations can only 
predict the leading logarithmic dependence of the distribution functions with Q 2 . These 
distributions have in addition uncalculable non-perturbative corrections suppressed by 
inverse powers of Q 2 , the so-called higher-twist contributions: 



F,(x,Q 2 )=F l LT (*,Q 2 ) + 



Q 2 



+ 



(6.30) 



The leading-twist term (LT) is the one predicted by perturbative QCD. Since the ad- 



37 



ditional l/{Q 2 ) n dependences have to be fitted from the data, they increase the final 
uncertainties. These corrections are numerically important for Q 2 < O(10 GeV 2 ) and for 
x close to 1. Obviously, the perturbative QCD predictions can be better tested at large 
Q 2 , where the higher-twist effects are smaller. 

Since the singlet structure functions are sensitive to the gluon distribution, which is 
badly known, they suffer from rather large errors. Good data at low values of x is needed 
in order to perform an accurate determination. The HERA experiments are making an 
important improvement in the knowledge of these distributions. The latest fits [ 53| , in- 
cluding the most recent HERA data, obtain gluon and sea-quark distributions at small 
x which are significantly different from those in previous standard sets of parton distri- 
butions. The new gluon distribution is larger for x < 0.01 and smaller for x ~ 0.1. The 
reduction of the gluon distribution in the interval x ~ 0.1 — 0.2 is compensated by an 
increase in the fitted value of a s [|53 , bringing the DIS determination [53|, 54 



a. (Ml) 



0.114 ±0.005 



(6.31) 



in better agreement with the world average values, which we discuss in the next section. 

7. DETERMINATION OF THE STRONG COUPLING 

In the massless quark limit, QCD has only one free parameter: the strong coupling 
a s . Thus, all strong interaction phenomena should be described in terms of this single 
input. The measurements of a s at different processes and at different mass scales provide 
then a crucial test of QCD: if QCD is the right theory of the strong interactions, all 
measured observables should lead to the same coupling. 
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Figure 28: Compilation of a s measurements as function of the energy scale |7], 
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Obviously, the test should be restricted to those processes where perturbative tech- 
niques are reliable. Moreover, the same definition of a s should be taken everywhere; the 
MS scheme is usually adopted as the standard convention. Since the running coupling is 
a function of energy, one can either compare the different determinations at the different 
scales where they are measured, checking in this way the predicted Q 2 dependence of the 
coupling, or use this prediction to bring all measurements to a common reference scale 
where they are compared. Nowadays, the Z-mass scale is conventionally chosen for such 
a comparison. 

In order to assess the significance of the test, it is very important to have a good 
understanding of the uncertainties associated with the different measurements. This is 
not an easy question, because small non-perturbative effects can be present in many ob- 
servables. In addition, some quantities have been computed to a very good perturbative 
accuracy (next-to-next-to- leading order), while others are only known at the leading or 
next-to-leading order; the resulting values of a s refer then to different perturbative ap- 
proximations. The estimate of theoretical uncertainties is also affected by the plausible 
asymptotic (i.e. not convergent) behaviour of the perturbative series in powers of a s . Al- 
though this is a common problem of Quantum Field Theories, it is probably more severe 
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in QCD because the coupling is rather big (at usual energies). 

Fig. 123 summarizes U\M the most reliable measurements of the strong coupling as 



function of the energy scale. The agreement with the predicted running of a s , indicated by 
the curves, is indeed very good. The value of a s (m%), extracted from the hadronic width 
of the t lepton, provides a very important low-energy measurement; although it has a 
rather large relative error, it implies a very precise prediction at the Mz scale, which is in 
excellent agreement with the direct determinations of a s (M|) performed at the Z peak. 
Fig. ^9] || compares the different measurements at the common reference scale Mz- The 



average af all determinations gives 

a s (Mf) = 0.117 ±0.005. (7.1) 



8. CHIRAL SYMMETRY 

Up to now, we have only discussed those aspects of QCD which can be analyzed in a 
perturbative way. Thus, we have restricted ourselves to the study of scattering processes 
at large momentum transfers, and inclusive transitions which avoid the hadronization 
problems. The rich variety of strong-interacting phenomena governed by the confinement 
regime of QCD has been completely ignored. 

There are certainly many approximate tools to investigate particular aspects of 
non-perturbative physics; however, rigorous first-principle QCD calculations seem unfor- 
tunately out of reach for present techniques. Nevertheless, we can still investigate some 
general properties of QCD using symmetry considerations. 

8.1 Flavour symmetries 

In order to build the QCD Lagrangian, we made extensive use of the SU(3)c colour 
symmetry, which is the basis of the strong interaction dynamics. The Lagrangian (|3.11|) 
has additional global symmetries associated with the quark flavour numbers: 
1- £qcd is invariant under a global phase redefinition of all quark flavours, 

Qf — * exp(i6)q f . (8.1) 

This symmetry is associated with the conservation of the baryon number. 

2. £qcd is also invariant under independent phase redefinitions of the different quark 
flavours, 

Qf — * exp(i6 f ) q f . (8.2) 

This symmetry implies the conservation of flavour. 

3. For equal quark masses, there is a larger symmetry under SU(Nf) transformations 
in flavour space, 

Qf — Uff Qf , Ue SU(N f ) . (8.3) 

This is a good symmetry of the light-flavour sector (u, d, s), where quark masses can 
be ignored in first approximation. One has then the well-known isospin (Nf = 2) and 
SU (3) symmetries. 

4. In the absence of quark masses, the QCD Lagrangian splits into two independent 
quark sectors, 

£qcd = ~\ G?G%, + iq L ^D^q L + iq R ^D^ R . (8.4) 



40 



Here, q denotes the flavour (and colour) vector q = column (it, d, .. .), and L, R 
stand for the left- and right-handed components of the quarks. Thus, the two quark 
chiralities live in separate flavour spaces which do not talk each other (gluon interac- 
tions do not change the chirality), implying that all previous flavour symmetries get 
duplicated in the two chiral sectors. 
The baryon number symmetry ( |8.1| ) is usually called Z7(l)y, since both chiralities 
transform in the same way. Its chiral replication is the corresponding U(1)a transforma- 
tion: 

q L — ► exp(-i9)q L ; q R — ► exp(i0) q R . (8.5) 

This symmetry of the classical (massless) QCD Lagrangian gets broken by quantum effects 
(triangular loops of the type shown in Fig. |5|, with gluons instead of photons); this is 
the so-called U(1)a anomaly. Although ( fg.5|) is not a true symmetry of QCD, it gets 
broken in a very specific way, which leads to important implications. A discussion of the 
phenomenological role of anomalies is beyond the scope of these lectures. However, let me 
mention that this anomaly is deeply related to interesting low-energy phenomena such as 
the understanding of the rj' mass, or the so-called proton spin crisi^. 

I want to concentrate here in the chiral extension of the old eightfold SU(3)v 
symmetry, i.e. in the global G = SU(3)l <8> SU(3)r symmetry of the QCD Lagrangian 
for massless u, d and s quarks. This larger symmetry is not directly seen in the hadronic 
spectrum. Although hadrons can be nicely classified in SU(3)v representations, degenerate 
multiplets with opposite parity do not exist. Moreover, the octet of pseudoscalar mesons 
(ir,K,r)) happens to be much lighter than all other hadronic states. 

There are two different ways in which a symmetry of the Lagrangian can be realized. 
In the usual one (Wigner-Weyl), the ground state (the vacuum) is also invariant. Then, all 



physical states can be classified in irreducible representations of the symmetry group |55 
Certainly, the hadronic spectrum does not look like that, in the case of the chiral group. 

There is a second (Nambu-Golstone), more sophisticated, way to realize a symme- 
try. In some cases, the vacuum is not symmetric. The hadronic spectrum corresponds to 
energy excitations over the physical vacuum and, therefore, will not manifest the original 



symmetry of the Lagrangian. However, Goldstone's theorem |56[ says that in such a case 
there should appear a massless scalar for each broken generator of the original symmetry 
group. If the chiral symmetry is realized in this way, there should be eight pseudoscalar 
massless states (Goldstone bosons) in the hadronic spectrum; this is precisely the number 
of states of the lightest hadronic multiplet: the 0~ octet. Thus, we can identify the it, 
K and rj with the Goldstone modes of QCD; their small masses being generated by the 
quark-mass matrix which explicitly breaks the global chiral symmetry of the Lagrangian. 

In the Standard electroweak model, the local SU (2)l<S>U(1)y symmetry is also real- 
ized in the Nambu-Goldstone way. There, the symmetry-breaking phenomena is assumed 
to be related to the existence of some scalar multiplet which gets a vacuum expectation 
value. Since a local symmetry gets (spontaneously) broken in that case, the Goldstone 



4) This is a quite unfortunate name, because: 1) the underlying QCD dynamics has little to do 
with the parton-model description of the proton spin; and 2) it is certainly not a crisis but 
rather a success of QCD. The failure of the naive quark-model description of an observable 
where gluons are predicted to play a crucial role (the anomaly), is indeed a clear experimental 
confirmation of the QCD dynamics. 
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modes combine with the gauge bosons giving massive spin-1 states plus the Higgs parti- 
cle. The QCD case is simpler, because it is a global symmetry the one which gets broken. 
However, something should play the role of the electroweak scalar field. Since quarks are 
the only fields carrying flavour, they should be responsible for the symmetry breaking. 
The simplest possibility is the appearance of a quark condensate 

v = (0\uu\0) = (0\dd\0) = (0|ss|0) < , (8.6) 

generated by the non-perturbative QCD dynamics. This would produce a dynamical 
breaking of chiral symmetry, keeping at the same time the observed SU(3)v symmetry. 

8.2 Effective Chiral Lagrangian 

The Goldstone nature of the pseudoscalar mesons implies strong constraints on 
their interactions, which can be most easily analyzed on the basis of an effective La- 
grangian. The Goldstone bosons correspond to the zero-energy excitations over the quark 
condensate; their fields can be collected in a 3 x 3 unitary matrix U((f>), 

(Okl^lO) — (8.7) 

which parametrizes those excitations. A convenient parametrization is given by 

f/(0)EEexp(V2$//), (8.8) 



where 



-4= + ^ K° 
V V6 / 



(8.9) 



The matrix U(<f>) transforms linearly under the chiral group, [gL,R £ SU($)l,r] 

Ql gLqL, qn gRqR => U{(j>) g R U{<f)) g[ , (8.10) 

but the induced transformation on the Goldstone fields <fi is highly non-linear. 

Since there is a mass gap separating the pseudoscalar octet from the rest of the 
hadronic spectrum, we can build a low-energy effective field theory containing only the 
Goldstone modes. We should write the most general Lagrangian involving the matrix 
U ((/>), which is consistent with chiral symmetry. Moreover, we can organize the Lagrangian 
in terms of increasing powers of momentum or, equivalently, in terms of an increasing 
number of derivatives (parity conservation requires an even number of derivatives): 

£es(U) = ]T C 2n . (8.11) 

n 

In the low-energy domain, the terms with a minimum number of derivatives will dominate. 

Due to the unitarity of the U matrix, UW = 1, at least two derivatives are required 
to generate a non-trivial interaction. To lowest order, the effective chiral Lagrangian is 
uniquely given by the term 

C 2 = ?-Tr[d ll tf& i u] . (8.12) 
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Expanding U(4>) in a power series in $, one obtains the Goldstone's kinetic terms 
plus a tower of interactions involving an increasing number of pseudoscalars. The require- 
ment that the kinetic terms are properly normalized fixes the global coefficient / 2 /4 in 
(|8.12| ) . All interactions among the Goldstones can then be predicted in terms of the single 
coupling /: 



C 2 = - Tr [dpQPQ] + ^ Tr 



($ <9 M $) (* &> $) 



(8.13) 



To compute the tttt scattering amplitude, for instance, is now a trivial perturbative 



exercise. One gets the well-known |57j Weinberg result [t = (p' + — p + 



Tt 



71 71 



71 71 



t/f 



(8.14) 



Similar results can be obtained for tttt — > Att, 6tt, 8tt, . . . The non-linearity of the effective 
Lagrangian relates amplitudes with different numbers of Goldstone bosons, allowing for 
absolute predictions in terms of /. Notice that the Goldstone interactions are proportional 
to their momenta (derivative couplings). Thus, in the zero-momentum limit, pions become 
free. In spite of confinement, QCD has a weakly- interacting regime at low energies, where 
a perturbative expansion in powers of momenta can be applied. 

It is straightforward to generalize the effective Lagrangian ( |8.12| ) to incorporate 
electromagnetic and semileptonic weak interactions. One learns then that / is in fact the 



pion-decay constant / ~ /„■ — 92.4 MeV, measured in n — > [iv^ decay |p8| . The corrections 
induced by the non-zero quark masses are taken into account through the term 



r 

*—"rr 



Tr 



M(U + U ] ) 



M = diag(m u ,m d ,m s 



.15) 



which breaks chiral symmetry in exactly the same way as the quark-mass term in the 
underlying QCD Lagrangian does. Eq. flS.15|) gives rise to a quadratic pseudoscalar-mass 
term plus additional interactions proportional to the quark masses. Expanding in powers 
of $ (and dropping an irrelevant constant), one has 



r 



P 



-Tr 



Qf 4 



:Tr 



C($ 6 // 6 ^ 



(8.16) 



The explicit evaluation of the trace in the quadratic mass term provides the relation 
between the physical meson masses and the quark masses: 



M 2 



K° 



M. 



(m u + m d ) — 



(m u + m d ) 



(m u + m s ) 



P 

\v\ 

I 2 



e + 0(e 2 



1.17) 



(m d + m s )^, 

l M 

~(m u + m d + 4m s )y + e + 0(e 2 ) 
3 J 
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where 

2f 2 (2m s -m u - m d ) 

Chiral symmetry relates the magnitude of the meson and quark masses to the size of 
the quark condensate. Taking out the common \v \/f 2 factor, Eqs. ( 8.17 ) imply the old 
Current Algebra mass ratios, 



Ml ± M 2 K+ M K o 3M 2 8 



(8.19) 



m u + m d (m u + m s ) (m d + m s ) (m u + m d + 4m s ) 

and [up to 0(m u — m d ) corrections] the Gell-Mann-Okubo mass relation 

3MJ = 4M| - Ml . (8.20) 

Although chiral symmetry alone cannot fix the absolute values of the quark masses, 
it gives information about quark-mass ratios. Neglecting the tiny 0{e) effects, one gets 
the relations 

m d -m u = (Mj <0 - M 2 K+ ) - (Ml Q - Ml + ) = Q ^ ^ 
m d + m u M%a 

2m r m "~ mrf = Mk^Mk = 12.6 . (8 . 22) 
2(m u + m d ) M\ 



In ( 8.21 ) we have subtracted the pion square-mass difference, to take into account the 
electromagnetic contribution to the pseudoscalar-meson self-energies; in the chiral limit 
(m u = m d — m s — 0), this contribution is proportional to the square of the meson 
charge and it is the same for K + and 7r + . The mass formulae ( |8.21| ) and ( |8.22| ) imply the 
quark-mass ratios advocated by Weinberg: 

m u :m d :m s = 0.55 : 1 : 20.3 . (8.23) 

Quark- mass corrections are therefore dominated by m s , which is large compared with m u 
and m d . Notice that the difference m d — m u is not small compared with the individual 
up- and down-quark masses; in spite of that, isospin turns out to be an extremely good 
symmetry, because isospin-breaking effects are governed by the small ratio (m d — m u )/m s . 

The $ 4 interactions in (|8.16 ) introduce mass corrections to the nir scattering am- 
plitude ( ETT3D , 

t-Ml 

T 2 

J 7T 



T ^ 



Since / ~ f n is fixed from pion decay, this result is now an absolute prediction of chiral 
symmetry. 

The lowest-order chiral Lagrangian encodes in a very compact way all the Current 
Algebra results obtained in the sixties |59|] . The nice feature of the chiral approach is its 



elegant simplicity, which allows to estimate higher-order corrections in a systematic way. 
A detailed summary of the chiral techniques and their phenomenological applications can 



be found in Ref. p8| . 
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9. SUMMARY 

Strong interactions are characterized by three basic properties: asymptotic freedom, 
confinement and dynamical chiral symmetry breaking. 

Owing to the gluonic self-interactions, the QCD coupling becomes smaller at short 
distances, leading indeed to an asymptotically-free quantum field theory. Perturbation 
theory can then be applied at large momentum transfers. The resulting predictions have 
achieved a remarkable success, explaining a wide range of phenomena in terms of a single 
coupling. The running of a s has been experimentally tested at different energy scales, 
confirming the predicted QCD behaviour. 

The growing of the running coupling at low-energies makes very plausible that 
the QCD dynamics generates the required confinement of quarks and gluons into colour- 
singlet hadronic states. A rigorous proof of this property is, however, still lacking. At 
present, the dynamical details of hadronization are completely unknown. 

Non-perturbative tools, such as QCD sum rules and lattice calculations, provide 
indirect evidence that QCD also implies the proper pattern of chiral symmetry breaking. 
The results obtained so far support the existence of a non-zero q-q condensate in the QCD 
vacuum, which dynamically breaks the chiral symmetry of the Lagrangian. However, a 
formal understanding of this phenomena has only been achieved in some approximate 
limits. 

Thus, we have at present an overwhelming experimental and theoretical evidence 
that the SU (3)c gauge theory correctly describes the hadronic world. This makes QCD the 
established theory of the strong interactions. Nevertheless, the non-perturbative nature 
of its low-energy limit is still challenging our theoretical capabilities. 
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APPENDIX A: SU(N) ALGEBRA 

SU(N) is the group of A x A unitary matrices, UW = WU = 1, with det U — 1. 
The generators of the SU (A) algebra, T a (a — 1, 2, . . . , A 2 — 1), are hermitian, traceless 
matrices satisfying the commutation relations 

[T a ,T b ] = if abc T c , (A.l) 

the structure constants f abc being real and totally antisymmetric. 

The fundamental representation T a = X a /2 is iV-dimensional. For N = 2, X a are the 
usual Pauli matrices, while for A — 3, they correspond to the eight Gell-Mann matrices: 



A x = 
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a73 



V 



They satisfy the anticommutation relation 



{X a ,X b } = j^S ab I N + 2d abc X c , 



1^ 



1 

(A.2) 

i o o N 
o 1 o 

0-2; 

(A.3) 



where In denotes the A-dimensional unit matrix and the constants d abc are totally sym- 
metric in the three indices. 

For SU(3), the only non-zero (up to permutations) f abc and d abc constants are 

1 



if 123 = f 147 = _f 156 — f 246 — f 257 — f 345 — _ f 367 - _L f 458 = _L_ f 678 _ 
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The adjoint representation of the SU (A) group is given by the (A 2 — 1) x (A 2 — 1) 
matrices (T A ) bc = —if abc . The relations 



Tr(A a A 6 ) =4T F 5 afe , 



Ti(T a A T b A ) = r d f bcd = C A 5 ab , 



C F = 



A 2 



2N 1 
C A = N, 

define the SU(N) invariants T F , C F and Ca- Other useful properties are: 

2 
A 

4 



(A.5) 
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Tr (A a A 6 A c ) = 2(d abc + if abc ) 
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APPENDIX B: GAUGE-FIXING AND GHOST FIELDS 

The fields G^ have 4 Lorentz degrees of freedom, while a massless spin-1 gluon has 
2 physical polarizations only. Although gauge invariance makes the additional degrees of 
freedom irrelevant, they give rise to some technical complications when quantizing the 
gauge fields. 

The canonical momentum associated with G%, = SCqqv/ S^OqG^) = G^ , 

vanishes identically for = 0. The standard commutation relation 

[G^x),Ul(y)}5(x°-y ) = ig^8^(x-y), (B.l) 

is then meaningless for /j, = v = 0. In fact, the field G° a is just a classical quantity, since 
it commutes with all the other fields. This is not surprising, since we know that there 
are 2 unphysical components of the gluon field, which should not be quantized. Thus, we 
could just impose two gauge conditions, such as G Q a = and VG a = 0, to eliminate the 2 
redundant degrees of freedom, and proceed working with the physical gluon polarizations 
only. However, this is a (Lorentz) non-covariant procedure, which leads to a very awkward 
formalism. Instead, one can impose a Lorentz-invariant gauge condition, such as d^G^ = 0. 
The simplest way to implement this is to add to the Lagrangian the gauge-fixing term 

£gf = ~(PGl) (d v G v a ) (B.2) 

where £ is the so-called gauge parameter. The 4 Lorentz components of the canonical 
momentum 

= = G% - i g, {d u G a u ) (B.3) 



are then non-zero, and one can develop a covariant quantization formalism. Since (|B.2|) 
is a quadratic G% term, it modifies the gluon propagator: 

r H 4 k p- ik ( x -y) ( h^h v 1 

(oirKWG^no)^^/^^-!-^ + (1-0^1 ■ (B.4) 

Notice, that the propagator is not defined for £ = oo, i.e. in the absence of the gauge-fixing 



term (|R2|) 



In QED, this gauge- fixing procedure is enough for making a consistent quantization 
of the theory. The initial gauge symmetry of the Lagrangian guarantees that the redun- 
dant photon polarizations do not generate non-physical contributions to the scattering 
amplitudes, and the final results are independent of the arbitrary gauge parameter £. In 
non-abelian gauge theories, like QCD, a second problem still remains. 

Let us consider the scattering process qq — > GG, which proceeds through the 



three Feynman graphs shown in Fig. 30. The scattering amplitude has the general form 



T = J^'e^e^ . The probability associated with the scattering process 

v ~ E4 A) 4 A) * (b-5) 

Z A A' 

involves a sum over the final gluon polarizations. One can easily check that the physical 
probability Vt, where only the two transverse gluon polarizations are considered in the 
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q G 




q G 



a) b) c) 

Figure 30: Tree-level Feynman diagrams contributing to qq — > GG. 



sum, is different from the covariant quantity Vc, which includes a sum over all polarization 
components: Vc > Vt- In principle, this is not a problem because only Vt has physical 
meaning; we should just sum over the physical transverse polarizations to get the right 
answer. However, the problem comes back at higher orders. 




Figure 31: 1-loop diagrams contributing to qq — > qq. 



The covariant gluon propagator ( |B.4j ) contains the 4 polarization components; 
therefore higher-order graphs such as the ones in Fig. |31] get unphysical contributions 
from the longitudinal and scalar gluon polarizations propagating along the internal gluon 
lines. The absorptive part of these 1-loop graphs (i.e. the imaginary part obtained putting 
on-shell the two gluon lines within the loop) is equal to \T(qq — > GG)\ 2 . Thus, these loops 
suffer the same probability problem than the tree-level qq — > GG amplitude. The propa- 
gation of unphysical gluon components implies then a violation of unitarity (the two fake 
polarizations contribute a positive probability). 

In QED this problem does not appear because the gauge-fixing condition d^A^ = 
still leaves a residual gauge invariance under transformations satisfying DO = 0. This 
guarantees that (even after adding the gauge-fixing term) the electromagnetic current is 
conserved, i.e. <9 M J^ m = d^eQ^fj^) = 0. If one considers the e + e~ — > 77 process, which 
proceeds through diagrams identical to a) and b) in Fig. |3(], current conservation implies 
kpJ^ = AyJ w *' = 0, where fc M and hf , are the momenta of the photons with polarizations 
A and A', respectively (remember that the interacting vertices contained in jw' are in fact 
the corresponding electromagnetic currents). As a consequence, the contributions from the 
scalar and longitudinal photon polarizations vanish and, therefore, Vc = Vt- 

The reason why Vc 7^ Vt in QCD stems from the third diagram in Fig. involving 
a gluon self-interaction. Owing to the non-abelian character of the SU(3) group, the gauge- 
fixing condition d^G^ = does not leave any residual invariance Q. Thus, fc^jw*' ^ 0. 

5) To maintain d^{G^)' = after the gauge transformation (l3.8|) , one would need O50 a = 
g s f abc Q^(§0 b s )Q^ which is not possible because is a quantum field. 
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m x ♦ q x $ x q 

d) d-) 

Figure 32: Feynman diagrams involving the ghosts. 

Again, the problem could be solved adopting a non-covariant quantization where 
only the physical transverse polarizations propagate; but the resulting formalism would be 
awful and very inconvenient for performing practical calculations. A more clever solution 
consist |6(J in adding additional unphysical fields, the so-called ghosts, with a coupling 
to the gluons such that exactly cancels all unphysical contributions from the scalar and 
longitudinal gluon polarizations. Since a positive fake probability has to be cancelled, 
one needs fields obeying the wrong statistics (i.e. of negative norm) and thus giving 
negative probabilities. The magic cancellation is achieved by adding to the Lagrangian 



the Faddeev-Popov term |6T 



C FP = -d^D^r , D^ a = d»r - g s f abc b G» , (B.6) 

where a , a (a = 1, . . . , N% - 1) is a set of anticommuting (i.e. obeying the Fermi-Dirac 
statistics), massless, hermitian, scalar fields. The covariant derivative D^<fi a contains the 
needed coupling to the gluon field. One can easily check that diagrams d) and d') in Fig. [32] 
exactly cancel the unphysical contributions from diagrams c) and c') of Figs. [30] and [31], 
respectively; so that finally Vc — "Pt- Notice, that the Lagrangian (P-6Q is necessarily not 
Hermitian, because one needs to introduce an explicit violation of unitarity to cancel the 
unphysical probabilities and restore the unitarity of the final scattering amplitudes. 

The exact mechanism giving rise to the £fp term can only be understood (in a 
simple way) using the more powerful path-integral formalism, which is beyond the scope 
of these lectures. The only point I would like to emphasize here, is that the addition of the 
gauge-fixing and Faddeev-Popov Lagrangians is just a mathematical trick, which allows 
to develop a simple covariant formalism, and therefore a set of simple Feynman rules, 
making easier to perform explicit calculations. 
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